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Abstract 

The role of f3 decay as a low-energy probe of physics beyond the Standard 
Model is reviewed. Traditional searches for deviations from the Standard 
Model structure of the weak interaction in /3 decay are discussed in the 
light of constraints from the LHC and the neutrino mass. Limits on the 
violation of time-reversal symmetry in (3 decay are compared to the strong 
constraints from electric dipole moments. Novel searches for Lorentz sym¬ 
metry breaking in the weak interaction in [3 decay are also included, where 
we discuss the unique sensitivity of /3 decay to test Lorentz invariance. We 
end with a roadmap for future /3-decay experiments. 
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I. INTRODUCTION 


The study of nuclear and neutron (3 decay has played a major role in uncovering the structure 
of the weak interaction, and therefore in the development of the electroweak sector of the 
Standard Model (SM) of particle physics. The intensity and the variety of {3 emitters, 
combined with the high precision with which /S-decay parameters can be measured, ensured 
that (3 decay remained important in searches for new physics beyond the SM (BSM). Novel 


techniques of laser cooling and atom trapping (Behr and Gwinner, 2009 Sprouse and Orozco 


1997) made it possible to detect the momentum of the recoiling nucleus, allowing for further 


searches in unexplored observables that became available. New sources for slow neutrons 


enabled further progress in the study of neutron /1-decay observables (Abele, 2008 Dubbers 


and Schmidt, 2011 Nico, 2009). The motivation for these modern experiments is on the 


one hand to improve the accuracy of SM parameters, and on the other hand to search for 
physics BSM. 

Searches for BSM physics in f3 decay look for deviations from the left-handed vector- 


axial-vector (“U — A”) space-time structure of the weak interaction, see Holstein| (2014) 


and Severijns et al. (2006) and references therein. High-precision /1-decay experiments are 
sensitive to possible contributions of non-SM (or exotic) currents, in particular right-handed 
vector, scalar, and tensor currents, that couple to hypothetical new, heavy particles. These 
exotic currents can also give additional violations of the discrete symmetries parity (P), 
charge conjugation (C), and time-reversal invariance (T). 

Traditionally, (3 decay has been viewed as complementary to the direct searches for new, 
heavy particles at high-energy colliders. However, with the availability of meson factories 
the emphasis of searching for new physics in precise measurements of semileptonic decay 
parameters has shifted from (3 decay. New physics has also been severely constrained by the 
emergence of the new field of neutrino oscillations and by the ultra-precise measurements 
of static observables such as the weak charges of quarks and electrons and the P- and T- 
odd electric dipole moments (EDM) of particles, atoms, or molecules. Moreover, theoretical 
developments made it clear how various observables are interconnected, and therefore how 
the discovery potential of /9-decay experiments compares to that of other helds. 

Recently, another twist has been added to (3 decay as a promising precision laboratory to 
test the invariance of the weak interaction under Lorentz transformations, that is, boosts and 
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rotations. The available evidence for the Lorentz invariance of the weak interaction is, in fact, 
surprisingly poor. The possibility to break Lorentz and the closely related CPT invariance 


(Greenberg, 2002) occurs in many proposals that attempt to unify the SM with general 
relativity, one of the central open issues in theoretical high-energy physics. During the last 
decade, the phenomenological consequences of such a breakdown of Lorentz symmetry have 


been charted (Colladay and Kostelecky, 1998), and recently such theoretical studies have 


been extended to (3 decay (Noordmans et al. ^ 2013b). 

This review gives a broad overview of the searches for symmetry violations in nuclear 
and neutron [3 decay and discusses their signihcance compared to various other observables, 
both in precision measurements and in collider searches. In this way, it attempts to identify 
which /3-decay studies are the most relevant to pursue. In Sec. |TT] we hrst introduce the 
effective held theory (EFT) framework, which enables us to compare various experiments in 


a model-independent approach. We dehne the /3-decay observables in Sec. Ill 


In Sec. |IV| we review the best bounds on exotic right-handed vector, scalar, and tensor 
couplings. We hrst address the most sensitive /3-decay experiments, in which we also include 
limits from pion-decay experiments. 

Second, we discuss how the neutrino mass and data from the LHC experiments constrain 
BSM physics. We compare the bounds from these two sectors with the bounds from /3-decay 
experiments. The violation of time-reversal invariance is discussed in Sec. |Vj In /3 decay, 
T-violation manifests itself in nonzero imaginary parts of the couplings, which are probed 
by triple-correlation observables in [3 decay. We discuss how these bounds compare to those 
derived from the stringent upper limits on the values of EDMs. 


In Sec. VI, we address the possibility that the weak interaction violates Lorentz symme¬ 
try, and in particular rotational invariance, in nuclear and neutron (3 decay. Such Lorentz 
violation (LV) would give rise to unique signals with no SM “background,” which, even when 
tiny, could be experimentally detectable. Nuclear and neutron (3 decay offer a unique sen¬ 
sitivity to some Lorentz-violating parameters, especially in the gauge and neutrino sector, 
which we discuss separately. 

We conclude with a roadmap for the opportunities in future /3-decay studies, in the light 
of the obtained and foreseen bounds from other frontiers. 
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Table I Classification of nuclear (5 decays and their characteristic use in the SM and in the search 
for BSM physics. 


II. FORMALISM 


Nuclear and neutron (3 decay are semileptonic processes, mediated by the W gauge boson 
of the electroweak interaction. This interaction is described by a spontaneously broken 
SU{2)l X U{1)y gauge symmetry. Under SU{2)l symmetry, left-handed leptons transform 
as a doublet, while right-handed particles are SU{2)i singlets. This is denoted by 

La = ) LIa = {Ia)r , (1) 

where A is the flavor index and the left- and right-handed fields are 

V’L = ^( 1 -75)^ , V’i? = ^( 1 + 75)^^ • ( 2 ) 

The W boson only interacts with left-handed fermions, which reflects the maximal violation 
of parity (P) symmetry in the weak interaction. In the minimal SM neutrinos are assumed 
to be massless, and right-handed neutrinos are absent. The role of the neutrino mass is 
discussed in Sec. IIV.CI 

The decay transition d —)■ ue~i)e {u —)■ de'^Ue) is, in the limit of inhnite lU-boson 

mass, described by the effective Lagrange density 

- 75)«^e h7^(l - lb)d + h.c. , (3) 
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where Gp is the Fermi coupling constant, Vud is the ud entry of the Cabibbo-Kobayashi- 
Maskawa (CKM) mixing matrix, and h.c. denotes the hermitian conjugate. We work in 
natural units, h = c = 1, and use 7 ® = and = —60123 = 1- 

At the nucleon level, all possible quark bilinears and their associated form factors need 


to be inserted (Weinberg, 1958), such that 


(p|My^d|n) = p 

{p\uiiil5d\n) = p 


9v{q^)li. + 


9M{q^) 


9A{q^)l^,'lb + 


M 

9T{q^ 

2M 


(^t,uq‘' + 


~9s{q^ 


2M 


-qf^ 


-(T^,uq''lb + 


~9p{q^ 

2M 


n 


-g^75 


n , 


( 4 ) 


where q = Pn — Pp is the momentum transfer and M is the nucleon mass. The vector form 
factor gv and the axial-vector form factor qa give the leading contributions to f3 decay, 
because the nuclei can be treated nonrelativistically. In the isospin limit, the induced form 
factor qm, called weak magnetism, is given by {pp — Pn)/‘^, he. the difference between the 
magnetic moments of the proton and the neutron. Given the current experimental precision, 
this form factor can be neglected, but future experiments might reach a level of precision 


for which weak magnetism has to be taken into account, see Sec. IV.D In the isospin limit 


the induced scalar form factor gs and tensor form factor gp vanish (Weinberg, 1958), and 


we can neglect them at present. The induced pseudoscalar form factor gp gets an additional 
suppression of q/M, because of the pf^n structure. We comment on pseudoscalar couplings 
in Sec. IIV.A.5I 

The leading-order SM expression for neutron decay is 


/IsM = gviq^)e'yp{l - y 5 )i/ePy''(l - + h.c. 


V2 


9v 


( 5 ) 


In the limit of —)■ 0, the vector charge is (7v(0) = 1, up to small corrections. This is 
dictated by the hypothesis of the conserved vector current (CVC). The axial-vector charge 
gA is only partially conserved (PCAC). The best current value is derived from neutron 


/3-decay experiments, \gA\ = 1.2723(23) (Olive et al.. 2014). 


In nuclear /3 decay one can exploit the properties of the parent and daughter nucleus 
to select particular parts of the interaction. Pure Fermi (F) transitions probe the vector 
currents (y^), while pure Gamow-Teller (GT) transitions probe the axial-vector currents 
( 757 ^)- Mixed transitions always require knowledge of the Fermi and Gamow-Teller tran- 
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sition matrix elements, Mp = (/|1|*) and Mgt = respectively. The conditions for 

spin change (A J) and parity change (vTjVrj) for Fermi and Gamow-Teller transitions are given 
in Table This Table also lists for which aspect in SM and BSM research these transitions 
are used. We have dehned the Fermi-Gamow-Teller mixing ration 

P = QaMgt/ dvMp , (6) 


and 


-^ = \9A\/gv 


(7) 


It is desirable to reduce the uncertainties of nuclear structure and select the simplest isotopes. 
For Fermi transitions the superallowed 0+ —)■ O’*" transitions are of most interest. For mixed 
transitions, mirror nuclei are preferred. For general mirror nuclei p has to be measured, 
while neutron decay (J^ = 1/2+ —)■ J'^ = 1/2+, \Mp\‘^ = 1 and \Mgt\^ = 3) allows for the 
determination of the value of A (Abele, 2008 Dubbers and Schmidt, 2011| Nico, 2009). An 


elaborate compilation of neutron-decay amplitudes is given in Ivanov et al. (2013). 


When searching for physics BSM, nuclei serve as “micro-laboratories” that can be judi¬ 
ciously chosen to look for certain manifestations of new physics. In this review, we address 
both the traditional searches for exotic couplings and the novel searches for Lorentz violation. 
In the latter, the possibility of angular-momentum violation needs to be considered, where 
the simplest of the forbidden decays, hrst-forbidden unique transitions, become relevant 


(Noordmans et al., 2013a). Both helds search for BSM physics generated by an unknown 


fundamental theory at a high-energy scale. To study the effect of new physics at low ener¬ 
gies, we work in an EFT approach. Within this framework the effects of new physics at low 
energies are described in a model-independent way with an effective Lagrangian of the form 


= -CsM + c 


-BSM 


( 8 ) 


The search for exotic couplings focuses on right-handed vector, scalar, and tensor couplings. 
These non-SM interactions can be included in the Lagrangian by adding higher-dimensional 
operators to /^bsm- The effects of Lorentz violation can also be described in an EFT frame¬ 


work (Golladay and Kostelecky, 1998 Noordmans et al, 2013b). We discuss both frame¬ 


works separately. 

























A. Exotic couplings 


In EFT, deviations from the V — A structure due to exotic couplings are generated by 
higher-dimensional operators, which are suppressed by the high-energy scale A. The effective 
Lagrangian is parametrized as 


= /:sM + , 


( 9 ) 


where 


£(4+t) = 


(4+fc) 


( 10 ) 


and where the c* are dimensionless constants and are dimension-(4 -|- k) operators. 

The SM only contains operators with mass dimension 3 or 4. For Lorentz-symmetric BSM 
physics, the lowest term we could add is There is, however, only one dimension-5 

operator, namely the operator that generates Majorana neutrino masses ( |Weinbei^ , 1979). 
In searches for exotic couplings we assume the neutrino mass to be small, and therefore 
we neglect this operator. We focus only on as even higher-dimensional terms are 

suppressed by additional powers of the large scale A. 


The that contribute to semileptonic charged decays are listed in Cirigliano et al 


(2013b, 2010). At low energies these dimension -6 operators generate the original vector 


(Cy), axial-vector (Ca), scalar (C 5 ), pseudoscalar (Cp), and tensor (Ct) couplings of Lee and 


Yang| (1956). At the quark level, the effective Lagrangian for {3 decay, with non-derivative 
four-fermion couplings, i^ 


£(eff) _ 


72 


\ ttes eYK ■ u'yf.ds 

e,S=L,R 


_CT 




+A^s evl ■ uds + «£ 

V2 V2 


( 11 ) 


where we sum over the chirality (L, R) of the hnal states. 

The coefficients represent 

• a^ 5 \ all possible V and A couplings, 

• A^s'- exotic scalar/pseudoscalar couplings (where e denotes the chirality of the neutrino 


We follow Herczeg (20011, except for a factor GplTd/v^ that we have extracted. 


9 





















and 6 the chirality of the d quark), 


a^: exotic tensor couplings (where e denotes the chirality of both the neutrino and the 
d quark). 


These coefficients are related to the couplings Ci and C' {i = S, V, A, T, P) of Lee and Yang 
(1956) by Eqs. (A3) and (A4) of Appendix [A} In the SM all couplings except ull = 1 are 
zero. For tensor couplings, only aL and aR occur, since cr^ul 5 = ■ The constants 

0 ^ 5 , A^s, and can be related to Cj, by matching their values at the low-energy scale with 
standard EFT techniques. The chiral structure of the coefficients is expressed by the first 
and second index, which denote the chirality of the neutrino and the d-quark, respectively. 
All couplings with first index R involve a right-handed neutrino. In the SM, right-handed 
neutrinos are absent, but they are present in many new-physics models. The role of the 
right-handed neutrino is discussed in Sec. |IV.C[ The new exotic couplings can be complex, 
representing the possibility of time-reversal (T) violation (Sec. 0. The introduction of left- 
handed and right-handed couplings leads to parity violation when the coefficients differ. In 
the absence of right-handed couplings, parity violation is maximal. 


To describe (3 decay of the nucleon we define the hadronic matrix elements (Herczeg 


2001 ) 


{p\u-f^d\n) = gv{q^)Pl,,n , 

(12a) 

{p\ul^,l5d\n) = gA{q'^)ini^'y5n , 

(12b) 

{p\ud\n) = gs{q‘^)pn , 

(12c) 

(p|M75d|n) = gp{q^)pi^n , 

(12d) 

{p\ua^ud\n) = gT{q^)p(Jf,un , 

(12e) 

Lgian in Eq. ([IT|) accordingly. As before, the vector charge is 


gy = gyiQ) = 1. The other couplings, gA,gs,gpy and gx can be calculated theoretically by 
using lattice QCD. Estimates for g^ on the lattice are currently not competitive with the ex¬ 


perimental value |5'yi(0)| = 1.2723(23) determined from neutron /3 decay (Olive et al., 2014). 
The scalar, pseudoscalar, and tensor constants, gs,gp, and gr, are determined theoretically. 


They are further discussed in Sec. IV 


Searches for exotic coupling also include searches for right-handed V + A currents. Such 
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currents are predicted for instance by left-right (LR) models, which add an SU{2)r gauge 
symmetry to the SM. This extends the SM with an additional gauge boson Wr, which 
mixes with the original SM W boson Wr. The weak eigenstates can be expressed in the 
mass eigenstates Wi and W 2 as 

Wl = Wi cos ^ -|- W 2 sin ^ , (13a) 

Wr = sin^ + W 2 cos^ , (13b) 


where ^ is the mixing angle and u a CP-violating phase. The coupling of Wr to quarks and 
leptons introduces the right-handed coupling qr and the right-handed CKM element the 
equivalents of the SM parameters. The expressions for aLR,aRL, and aRR in terms of these 


parameters are given in Herczeg (2001). A specihc class of LR models are the symmetric 
LR models, in which P or C symmetry of the Lagrangian is imposed, which implies Ql = Qr- 
We focus on bounds for such models in Sec. IIV.BI 


B. Lorentz violation 


The study of Lorentz violation is motivated by the possibility of spontaneous breaking 


of Lorentz invariance predicted by theories of quantum gravity (Kostelecky and Samuel 


1989 Liberati, 2013 Liberati and Maccione, 2009). The natural energy scale for these 


theories of quantum gravity is the Planck scale, which lies 17 orders of magnitude higher 
than the electroweak scale. This precludes the direct detection of Planck-scale physics, but 
the effects of Lorentz violation at the Planck scale can become manifest at much lower 
energies, providing a “window on quantum gravity.” At low energy, Lorentz violation can be 
systematically described by the Standard Model Extension (SME) ([Colladay and Kostelecky 


1998), by using an EFT approach. The SME contains all possible Lorentz-violating terms 


that obey the SM gauge symmetries, which include CPT-violating terms, since Lorentz 
violation allows for the breaking of CPT invariance. In fact, CPT violation can only occur 


if Lorentz symmetry is also broken (Greenberg, 2002). 


Spontaneous Lorentz violation arises as Lorentz-tensor helds acquire a vacuum-expectation 
value (VEV), resulting in Lorentz-violating tensor coefficients in the SME Lagrangian. These 
coefficients can be understood as constant background tensor helds. Due to these tensor 
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fields, the Lagrangian is no longer invariant under particle or active Lorentz transformations, 
i.e. boosts or rotations of the particles, because the background helds do not transform 


under the Lorentz group (Colladay and Kostelecky, 1998). However, the low-energy theory 
remains invariant under observer Lorentz transformations, i.e. boosts or rotations of the ob¬ 
server’s inertial frame. Because Lorentz symmetry is spontaneously broken, the underlying 
fundamental theory at the Planck scale remains Lorentz invariant, implying that important 
features such as energy-momentum conservation and microcausality are still valid. A possi¬ 
ble experimental signature of Lorentz violation is a sidereal variation of observables, which 
arise as the laboratory moves through the Lorentz-violating background held when Earth 


rotates (other examples are given in e.g. Mattingly (2005)). 

Schematically, terms in /Ibsm in Eq. (| 8 |) can be written as (Colladay and Kostelecky 

y(4+fc) 


1997) 


£np = (T) + — (T) ■ (T) ■ C>W+'=) . 


(14) 


where we summed over repeated indices and where are dimensionless constants, (T) is the 
expectation value of tensor T, P = 1, 75 , 7 ^, 7 ;_j 75 , represents the gamma-matrix structure, 
and are higher-dimensional operators. Furthermore, A represent the scale of the 

fundamental theory, which would naturally be the Planck scale. The higher-dimensional 


operators are suppressed by powers of this high scale. The hrst two terms in Eq. (14) have 
mass dimension 3 and 4, respectively. These terms are described in the original SME papers 


by Colladay and Kostelecky (1998) and are now referred to as the minimal Standard-Model 
Extension (mSME). For our present discussion we limit ourselves to the mSME, although 
higher-dimensional coefficients have also been described (Bolokhov and Pospelov[ (2008); 


Kostelecky and Mewes (2009, 2012, 2013)). 


From an EFT point of view, the introduced Lorentz-violating dimension-3 and dimension- 
4 operators are unnatural. Naively, one would expect the dimension-3 operators to scale 
linearly with the large scale A, while the coefficients of the dimension-4 operators should be 
of order unity. The experimental bounds on these dimension-3 and dimension-4 operators 
are much smaller, of course. This problem does not occur for higher-dimensional opera¬ 
tors, which are naturally suppressed by the scale A. To evade these naturalness problems, 
the current limits on dimension-3 and -4 coefficients require either large hne-tuning, or a 
symmetry that forbids these coefficients. However, even if dimension-3 and -4 operators are 
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forbidden at tree level, they will be induced by quantum corrections generated by higher¬ 
dimensional non-renormalizable operators. These corrections scale quadratically with the 
cutoff scale, which might be as large as A. This can be circumvented by introducing new 
physics between the weak scale and the Planck scale. In that case, radiative corrections 


scale with a signihcantly lower cutoff scale (see e.g Mattingly (2008)). Such a scenario oc¬ 


curs in supersymmetry (SUSY) (Bolokhov et al ., 2005; Groot Nibbelink and Pospelov, 2005). 
SUSY restricts Lorentz-violating operators to dimension 5 and higher, and forbids those of 
dimension 3 and 4. Dimension-3 and dimension-4 operators are generated by loop correc¬ 
tions if SUSY is broken. This would naturally lead to a suppression of m?/K and m/A for 
dimension-3 and dimension-4 operators, respectively, where m is the SUSY-breaking scale 


(Bolokhov et al. ^ 2005; Groot Nibbelink and Pospelov, 2005). In the mSME, it is assumed 
that dimension-3 and dimension-4 operators are suppressed by some unspecihed higher- 
scale mechanism, and the experimental constraints are studied without any assumptions on 


the nature of this suppression mechanism (Golladay and Kostelecky, 1998; Kostelecky and 


Russell, 2011). 


The SME contains a large number of coefficients that parametrize possible Lorentz vio¬ 
lation. We list the relevant coefficients for [3 decay, which are the lepton, Higgs, and gauge 


terms. The Lorentz-violating terms for leptons are (Golladay and Kostelecky, 1998) 


lepton = La i{cY^)L)^ “ 


^B 


+Ra i{c^)f,uAB'y''L)^ - {a'^)i,ABT 


LV\ 


R 


B 


(15) 


where L denotes the SU{2)l doublet and R denotes the singlet, defined in Eq. ([^. The 
subscripts A, B are flavor indices, and is the covariant derivative. This introduces the 
Lorentz-violating coefficients and which are GPT-odd and GPT-even, respectively. 
We have introduced the superscript LV for these coefficients, in order not to confuse them 


with the coefficients in Eq. (11). 


Before electroweak symmetry breaking, the Higgs and gauge sector are described by 
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(Colladay and Kostelecky, 1998) 


C 


Higgs+gauge 


+ H.c. + \i + H.c. 




2 >B 


fll/ 






( 16 ) 


where W^'', and are the SU{3)c, SU{2 )l, and U{1)y held-strength tensors, respec¬ 
tively, and 0 is the Higgs doublet. The coefficient is CPT-odd, and the only coefficient 
with dimension of mass. The other coefficients are CPT-even and dimensionless. The coef- 
hcient has symmetric real and antisymmetric imaginary components. The k^w and k^B 
coefficients are real and antisymmetric. The gauge couplings kc, kw, and ks are real and 


have the symmetry properties of the Riemann tensor (Colladay and Kostelecky, 1998). 


The SME parameters have been studied in a wide range of experiments (Kostelecky 


and Russell, 2011). The electromagnetic and gravity sector have been studied extensively. 


whereas the number of searches in the weak interaction is rather low. This changed recently 
(Miiller et al. ^ 2013 Noordmans et al., 2013a|[b ), and the search for Lorentz violation has 
been extended to weak decays, in particular (3 decay. {3 decay places strong constraints 
on Lorentz-violating coefficients in the Higgs and gauge sector. In addition, (3 decay has a 
unique sensitivity to some coefficients in the neutrino sector ( Diaz et 'ol] 2013). We discuss 
these constraints in Sec. lYD 


III. OBSERVABLES IN 3 DECAY 

A. Correlation coefficients in 3 decay 

In 3 decay, the correlations between different observables, such as the 3 momentum and 
the nuclear spin, can be measured. The amount of correlation is expressed in terms of 
correlation coefficients. These correlation coefficients depend on SM couplings and possible 
new V, A, S, P, and T interactions. Using the general effective Lagrangian in Eq. (0 , we 
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can write the decay-rate distribution for polarized nuclei as (Jackson et al.\ 1957b) 


u{{J)\Ee, fie, flu)dEedfledflu 


[2n) 


X \ 1 -|- (3 


Pe-Pv 

EeEy 


, rrie 


^ Pe-Pv 

3 EeEy 


(A 

J 


A^ + B^ + D 


E. 


E„ 


EeEi, 


{Pe-j){Pu-j) 

EeEy 


j(j + 1)-3((j.J)2) 

J(2J- 1) 


(17) 


where Ee(u), fle(u), and pe{u) denote the total energy, direction, and momentum, respec¬ 
tively, Eq is the energy available to the electron and the neutrino, (J) is the expectation value 
of the spin of the initial nuclear state, and j is the unit vector in this direction; E{±Z, Ee) 
is the Fermi function which modifies the phase space of the electron due to the Coulomb 
held of the nucleus. Also affecting the phase space is the Fierz interference term, factorized 
with the coefficient b. This term is zero in the SM. We dehned ^ where ^ gives 

the strength of the interaction. The remaining terms describe the /^-correlation coefficients: 
the /9-neutrino asymmetry a, the P-odd “Wu-parameter,” the /9-asymmetry A, the neutrino 
asymmetry B, and the triple-correlation coefficient D. The c coefficient vanishes for non- 
oriented nuclei and for nuclei with J = 1/2, such as the neutron. The c coefficient has not 
been taken into account in any experiment to date. However, in future experiments, which 
use laser beams to trap and cool samples, the expectation value ((J ■ j)^) may be affected, 
such that the c coefficient can play a role. 


The decay rate integrated over neutrino direction, but taking into account electron po¬ 


larization, is (Jackson et ai, 1957b) 


cn((t/), cTg IF/g, 


^^AA^ PeEeiEo - E.fdEedVle 

(27r)4 


J-CTg 


J ^ Ee + m \jEe J Ee 


, ( 18 ) 


where ag is the spin vector of the jd particle. This introduces the longitudinal jd polarization 
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Goefficient 

Gorrelation 

P 

T 

a (/9z/ angular correlation) 

Pe ' Pv / 

Even 

Even 

b (Fierz interference term) 

rrie/Ee 

Even 

Even 

A {P asymmetry) 

J ' Pe/E/e 

Odd 

Even 

B [u asymmetry) 

J -Pu/E^ 

Odd 

Even 

G (Longitudinal polarization) 

de ■ Pe/ Ee 

Odd 

Even 

N 

J ■ (Te 

Even 

Even 

Q 

de -PeJ ■ Ve/Ee 

Even 

Even 

D (triple correlation) 

J ■ (Pe ^ Pu) / EeE^, 

Even 

Odd 

R (triple correlation) 

ae-{J X Pe)/Ee 

Odd 

Odd 


Table II Overview of symmetry properties under parity (P) transformations and time reversal (T) 
of the most relevant correlation coefficients in allowed jS decay. 


G, the spin-correlation coefficients N and Q, and the triple-correlation coefficient R. The 
symmetry properties of the correlation coefficients are listed in Table [TT| The A, B, and 
G coefficients are associated with parity violation. Depending on the type of transition 
they can have SM values close to ±1, which is characteristic for maximal parity violation. 
The triple-correlation coefficients D and R are T-odd and unmeasurably small in the SM 


(Herczeg and Khriplovich, 1997). 


Integrating the decay rate over all kinematical variables gives the inverse lifetime, 




(19) 


where / contains the integration over the modihed phase space and {me/Ee) is the average 
inverse energy in units of the electron mass. 

In Appendix|^we list the relevant correlation coefficients in terms of the couplings dehned 


in Eq. (11) and the Fermi/Gamow-Teller matrix elements. The different correlation coeffi¬ 
cients contain combinations of the complex V, A, S, P, and T couplings. Given the current 
experimental precision, we have neglected Goulomb corrections. These corrections mainly 


introduce additional imaginary couplings (except for the D and R coefficients) (Jackson 


et al ., 1957b) 


We proceed by discussing how /S-decay correlation experiments, combined with lifetime 
measurements, are used to obtain precise values for the SM V and A coupling strengths. In 


Sec. IV we discuss constraints on exotic couplings. 
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B. Standard Model parameters in /3 decay 


The correlation coefficients in Appendix reduce to the SM expressions when putting 
the scalar and tensor couplings to zero, All^lr,rr,rl = 0 and aL(R) = 0, and by using only 
V — A couplings, aLFi,RR,RL = 0. The Fierz-interference coefficient b is zero in the SM. The 


lifetime in Eq. (19) can be derived from the ft value, using the measured half-life t instead 


of r. In the SM, 


m" 


ft 27r^ lii(2) 


GlVl,gl\MA\l + \p\‘) ■ 


( 20 ) 


The SM value for Gp is obtained from muon decay (Webber et al ., 2011). It is important to 


note that if one considers non-SM contributions these may influence muon decay as well. In 
principle, qa is calculable using lattice QCD, but as mentioned before, current lattice calcu¬ 
lations are not as accurate as values derived from experiments and henceforth A = \gA\/Qv 
is considered a free parameter. In general, Mp and Mgt need to be derived from nuclear 
model calculations. For superallowed Fermi transitions p = 0 and Mp = a/2, in the isospin 


limit. Hardy and Towner (2009) analyzed all available superallowed Fermi transitions, and 


derived a value for the ud CKM matrix element. Since the ft values of superallowed tran¬ 
sitions should be equal, a large number of measurements could be combined, leading to the 
most precise value of Vud = 0.97425(22) ( [Hardy and Towner , 2009). In the analysis, de¬ 
tails of the isotope-dependent nuclear-structure corrections on the matrix element Mp (e.g. 
isospin breaking) and the phase-space modihcations are also considered. The superallowed 


transitions also give the best bound on the Fierz coefficient b in Eq. (19) by considering the 


energy dependence of the lifetime (Sec. IV.A.1). 


The parameters A and Vud can also be determined from /3-decay correlations in neutron 
decay and from the neutron lifetime (Abele, 2008 Dubbers and Schmidt, 2011 [Mco 2009 


Wietfeldt and Greene, 2011|. The best current values are A = 1.2723(23) (Olive et al. 


2014) and Vud = 0.9742(12) (Dubbers and Schmidt, 2011). The latter is more than five 
times less precise, see also Fig. 22 in Dubbers and Schmidt ( 2011| . The strong Gamow- 
Teller dependence of neutron decay and the precision of the neutron-decay parameters is 
such that neutron decay also plays an important role in searches for tensor currents, as we 
will discuss in Sec. IIV.A.3I 

Another class of nuclei for which the nuclear structure is relatively well known are the 
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mirror nuclei (Severijns et al, 2008). Like neutron decay, mirror decays are mixed Fermi- 
Gamow-Teller transitions. Extraction of Vud from lifetime measurements requires knowledge 
of the mixing parameter p, such that an additional measurement of at least one of the corre¬ 


lation coefficients is necessary. Naviliat-Cuncic and Severijns (2009) hnd Vud = 0.9719(17), 


using 5 available transitions. The important structure corrections to Eq. (20) for mirror nu¬ 


clei have been evaluated (Severijns et al., 2008), in analogy to the work of Hardy and Towner 


(2009) for superallowed Fermi decays. This new class of nuclei will broaden the spectrum 
of data and remove any possible bias in selecting only superallowed Fermi transitions in the 
determination of Vud- Measurements with this motivation were undertaken. For example. 


Shidling et al. (2014), Broussard et al. (2014), and Triambak et al. (2012) have measured 


the lifetime of two relevant mirror nuclei, ^®Ne and We will not review the status of this 
held here, but comment on their relevance in limiting left-handed tensor couplings via the 
Fierz-interference term in the next section. It demonstrates that the contribution of nuclear 
physics to high-precision SM data goes hand in hand with the searches for new physics in (3 
decay. 


IV. CONSTRAINTS ON EXOTIC COUPLINGS 


(3 decay played an important role in establishing the V — A structure of the SM, initially 
eliminating to a large extent the possible contributions of scalar and tensor interactions. 
Modern searches in nuclear (3 decay consider again scalar and tensor currents as possible 


very small deviations from the SM due to new physics (see e.g. Severijns et al. (2006) and 
Severijns and Naviliat-Cuncic (2011[)). 


The searches in f3 decay are part of a much wider search in subatomic physics for new 
physics. Comparison between different searches has become possible in an EFT framework 
by using the effective Lagrangian in Eq. (In}. At the quark level the relations between 
different observables are clean, but at the nucleon level they involve the nuclear form factors 
9 A,gs,gp, and qt- Accurate values for these parameters are necessary in order to compare 
different limits. Recently, signihcant progress on the accuracy of both gs and qt has been 
reported. First results for gp are also available. The most precise value for gp is calculated 


with lattice QCD. Two recent results are from Green et al. (2012), gp = 1.038(16), and 


Bhattacharya et al. (2014), gp = 1.047(61). 
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The calculation method used in these works gives a much larger uncertainty for gs- Esti¬ 


mates range from gs = 0.72(32) (Bhattacharya et al, 2014) to gs = 1.08(32) (Green et al 


2012). A value for gs can also be derived using the CVC relation and lattice calculations 


(Gonzalez-Alonso and Gamalich, 2014), 


175 ( 0 ) = -= 1 . 02 ( 11 ) , 

dm„ 


( 21 ) 


where both (Gonzalez-Alonso and Gamalich, 2014) and dirig = 


TTid — mu (Golangelo et a/., 2011) are obtained separately from lattice calculations. However, 


the determination of gs with Eq. (21) might underestimate the error, because correlations 


between the numerator and denominator are neglected. Such errors could be avoided by 


calculating the ratio in Eq. (21) directly on the lattice. Further efforts to reduce the error 


for gs directly on the lattice are being pursued (Bhattacharya et al. ^ 2012, 2014). 


The pseudoscalar constant gp can be calculated by using the PGAG relation. Gombined 


with lattice QGD results (Gonzalez-Alonso and Gamalich, 2014) one hnds 


gpiO) = -^gA = 349(9) , 


mn 


( 22 ) 


where M = (Mp -|- Mn)/2 is the average nucleon mass and rriq = {mu + md)/2 = 3.42(9) 
MeV is the average light-quark mass determined on the lattice (Golangelo et al.. [2011 ). 


According to the PDG, mg = 3.5^o;2 MeV (Beringer et al ., 2012), which gives a much larger 


error, gp = OdOl^g. Nevertheless, this shows that the pseudoscalar form factor is of order 
0(10^). In (3 decay, pseudoscalar terms are generally neglected, because they only occur as 
higher-order recoil corrections. This surpresses pseudoscalar interactions compared to scalar 
and tensor interactions. The large value of gp cancels this suppression to a large extent, 
and /9-decay experiments may be sensitive to pseudoscalar couplings after all. There are, 
however, already strong constraints on pseudoscalar couplings from pion decay, as we discuss 
in Sec. IIV.A.51 

In the remainder of this Section we comment on searches for exotic couplings in (3 decay 


(Sec. IV.A), but considering only real couplings. We compare these results with constraints 


from the Large Hadron Gollider (LHG) experiments (Sec. IV.B) and due to the nonzero 


mass of the neutrino (Sec. IV.G). Bounds on imaginary couplings are discussed separately 
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in Sec. |Vl 


A. Constraints from /3 decay 


In nuclear jS decays, exotic couplings are mainly searched for in either pure Fermi or 
pure Gamow-Teller decays. Pure Fermi transitions depend on vector and possibly scalar 
couplings, while pure Gamow-Teller transitions depend on axial-vector and possibly tensor 
couplings. The use of mixed transitions is necessary when searching for interference terms. 
Preferred are isotopes with a relatively simple nuclear structure, e.g. mirror nuclei, or the 
neutron. We discuss the constraints from Fermi, Gamow-Teller, and mixed decays separately, 
focusing on the best current experimental data. We discuss the constraints on scalar and 
tensor couplings, while assuming no additional vector or axial-vector interactions. For a 


£t of the data including these interactions we refer to Severijns et al. (2006), where also a 


review of the experimental techniques is given. We discuss V + A couplings in Sec. IV.B 


Most /3-correIation coefficients are measured by constructing asymmetry ratios. For ex¬ 
ample, the /3 asymmetry is measured from the quantity 


A 


measured — 


N{t) - Nil) 
iv(t) + iv(;) ’ 


(23) 


where N ('() and N (J,) are the decay rates derived from measuring f3 particles in a particular 
detector while the polarization, P, of the nucleus changes sign. The arrows indicate the 


direction of polarization. The rates V('(),iV(4,) correspond to the integration of Eq. (17) 
over all unobserved degrees of freedom, which removes the dependence on the neutrino 
direction. In the numerator only the P-odd term remains, while in the denominator the odd 
term drops out. However, the Fierz interference term remains in the sum iV('() -|- so 

that 


A 


measured 


Ln F{±Z, EMEq - E,)^A\P\{pjE,) cosOJEjn, 
kn n±Z, Pe)Pe(Eo - Pe)^(l + b/E,)dE,dn, 

A\P\{l3eCOS 9e) 

1 + Kt 


( 24 ) 
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This implies that actually not the coefficient A is measured, but 


A = 


A 


i+ds 


( 26 ) 


The inverse average energy is approximated by 


me \ F(±z, EMEq - E^fdE, 

E, / F{±Z, E,)p,{Eo - E^fEJE, ’ 


(26) 


which depends on the specific isotope and the experimental setup. In principle, the average 
energy could also depend on the angular distribution (6*e). This makes it preferable that the 
analysis of {rrie/Ef.) is done and published together with the observed correlation coefficients. 
At present, many of the values for (mg/T^e) are derived by using the /3-energy threshold Emin 


(Pattie et al.^ 2013 Severijns et al., 2006 Wauters et al. ^ 2014) 


For the measured quantity X, X = a,A,B,G, etc., Eq. (25) applies. Except for B and 


N, the numerator of Eq. (25) depends only on the square of the coupling constants, while 


b has a linear dependence on left-handed couplings. In such cases one is most sensitive 
to b, and the measurement of X provides in the first place a measurement of the Fierz 
coefficient b. Therefore, the exact value of the {me/Eg) will become increasingly important 
with increasing experimental precision. 


1. Nuclear scalar searches 


Throughout the discussion of limits on scalar and tensor couplings, we will assume conven¬ 
tional left-handed vector couplings for the V-A part, such that an = 1, and aLR,RL,RR = 0. 


These and the other couplings are defined in Eq. (11). The notation is chosen such that 


the difference between the left-handed and right-handed coupling of the neutrino is empha¬ 
sized, i.e. for the scalar couplings Al = All + Alr (left-handed neutrino coupling) and 
Ali = Arr + Arl (right-handed neutrino coupling). Further details on the notation and 
some relevant expressions can be found in Appendix 
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For pure Fermi transitions 


i = 2\MF\^gi.{l + 


9s^ 


Al + A^, 


^bp — ±4:'j\MF\^gvgsAL , 


(27) 

(28) 


from Eq. (A9) and Eq. (All), where bp is the Fermi part of the Fierz coefficient b, the upper 


(lower) sign is for decays and 7 = Vl — Z‘^a^, with Z the atomic number of the 

daughter nucleus and a the fine-structure constant. For the positron-emitting superallowed 
O’*" —)■ 0+ Fermi decays 


m: 


ftp 27r3 1 n( 2 ) 


Gl gl \Mp\ 


1+ P 

\9v, 


Al + A] 


- 27 ( ^ > 


E, 


9v 


(29) 


Hardy and Towner (2009) obtained an average of all ft values, Et, after the appropriate 


corrections for radiative and nuclear-structure effects. The current best value of Vud is 
derived from Et^ assuming no exotic couplings. Allowing for scalar terms one can exploit 


(Hardy and Towner, 2005) the different values of {me/Ef) to put a stringent limit on the 
Fermi Fierz-interference coefficient (Hardy and Towner| 2009), 


bp = -0.0022(26) = -2 


^Ap 

9 V ^ 


l + lf{Al + Al) 


~ -2 —Al 
9v 


(30) 


Although bp is not sensitive to right-handed scalar currents, the value of Et is sensitive to 
these. In fact, the bound on right-handed couplings is more than an order of magnitude 
larger than that of left-handed couplings, such that both contributions to the Et values are 


of the same order, as can be seen in Eq. (29). Therefore, in searches for BSM physics one 


may not assume Vud as given by the PDG when such a search concerns also right-handed 
scalar terms. In the correlation coefficients, the value of Vud mostly drops out, but in limits 
derived from measured lifetimes the actual value of Vud is required. 


Constraints on right-handed scalar couplings can be extracted from the /3-z/-correlation 

as the lower bound and 


^exp \ 


coefficient a dehned in Eq. (AlO). We dehne 5- = \asM 
5+ = |a+p — asM\ as the upper bound, where the experimental value, at 90% conhdence 
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level (C.L.), lies between and a+p. Limits from a then give 


^ j \Al + A%] < , 


( 31 ) 


which gives a circular bound in the Al,Ar plane. Thus, the bound on Al and Ar would 
be the same, 


g^A 

9v 


LiR) 


< by 


(32) 


In practice experiments normalize the correlation to the total number of counts, and the 
absolute normalization is not measured. This means that in fact a is measured, as discussed 


below Eq. (23). In this way the Fierz-interference term b enters. The bounds remain circular. 


but the bound on Al changes to 


-S_ 


2'y{melEe) gv 


<^A,< 


2'y {rriejE^) 


(33) 


for and with opposite signs for (3 . 

Figure [T] shows the bounds from the best current experiments. The superallowed Fermi 


decays only constrain left-handed couplings and give a narrow vertical band (Hardy and 


Towner, 2009). The right-handed coupling Ar is constrained only by the correlations. 


and depends on the square root of the experimental error 5_. The most sensitive (3-u corre¬ 


lation measurements are from (Gorelov et al ., 2005) and ^^Ar (Adelberger et al ., 1999). 
We also include the recent measurement of the mirror nucleus ^^Na (Vetter et al., 2008), a 


mixed transition, where we have put tensor contributions to zero. In an earlier review this 


was erroneously shown with a bound as in Eq. (32) (Severijns and Naviliat-Cuncic, 2011). 


We show it because it is the hrst mixed transition available with such competitive preci¬ 
sion. The best current bounds on real scalar couplings from pure Fermi decays are found by 


minimalizing the x^"distribution of the bp from Eq. (30) and the measurements of the P-u 


correlation in (Gorelov et al., 2005) and ^^Ar (Adelberger et al ., [1999[ ) (Table [nT|). At 
90% G.L., 


-0.1 X 10“^ < —Al < 0.3 X 10“^ , 
9v 

-6 X 10-2 < —Ar < 6 X 10-2 . 

9v 


(34a) 

(34b) 
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Figure 1 (Color online) Bounds on left- and right-handed scalar couplings (90% C.L. 
0 + 


The narrow 

O"*" band is from superallowed Fermi transitions Eq. (30) (Hardy and Towner 2009). The 


ring-shaped boundaries are derived from P-u correlation measurements in (Gorelov et al 


2005) 

and (T' 

Ldelberger et al. 

2iNa 

(Vetter et al. 

2008 

) is given, 


1999), cf. Eq. (33). Also the bound from the mirror nucleus 


For Al the bound comes from the strong limit on the Fierz-interference term. The limit 
on Aji is less strong. Improving the bound on right-handed scalar couplings substantially 
is a daunting task: exploiting the forward-backward symmetry in the P-u correlation would 
require collecting 10^"^ events to reach a bound < 10“^ on QbAr. 


2. Nuclear tensor searches 


The nuclear Gamow-Teller matrix element Mqt can only be evaluated in the context 
of a nuclear model, because the spin of a nucleus is an observable, but the orbital angular 
momentum of a valence nucleon is not. For this reason Mqt cannot be evaluated sufficiently 
robustly to put a bound on the left-handed tensor couplings from ft values, as was done for 
the scalar coupling by using the superallowed Fermi decays. However, the Fierz-interference 
term will enter most observables via the normalization requirement discussed previously, cf. 


Eq. (25). The /9-asymmetry coefficient A in Gamow-Teller decays is a good example of this. 
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where 


A = 


A, 


GT 


1 + 6 , 


GT 


I me 


— ±Aj/ 


j'j 


-l + 8^ai -4 
9a 


gr /me 

\gA\°“-''\E, 


( 36 ) 


from Eq. (A13). Thus in the absence of Coulomb corrections one Ends that A becomes 


independent of and therefore only limits on gl can be obtained from A. Dehning the 
experimental bounds of A — Asm as before gives 


-h_ 


^ 9t ^ 
< - ;OCL < 


6 , 


Is^l A-yim^/Ee) 


(36) 


To obtain a bound on aji one can exploit the P-u correlation a. The result is similar to the 
result for a in Fermi decay. For /3~ Gamow-Teller decay osm = “1/3 and the bounds are 


I9t I 36— 

-ttH < \ —— , 

9a V 8 

3h- Ot ^ 

< ^< 


3^4 


A'y{me/Ee) \gA\ A'y{me/Ee) 


(37) 


The limits on tensor interactions can be improved by combining scalar and tensor 
searches. In particular, the left-handed tensor couplings can be further constrained by 
using the measurements of the Fermi and Gamow-Teller-transition ratio of the longitudinal 
(3 polarization. These measurements were performed in the first place to study the manifest 
left-right symmetric model (Garnoy et al., 1991t [Wichers et al. 1987), see also Sec. IV.B 


was 


The ratio of longitudinal polarizations (P, see Appendix of the emitted positrons 
measured in the systems 26y^jm^30p (Wichers et a/. , |1987 ) and (Garnoy et al, 1991), 

where the hrst nucleus decays via a Fermi and the second a Gamow-Teller transition. The 
two transitions have nearly identical endpoint energies, which eliminates systematic errors. 
The measured ratio is 


Pf 


P 




GT 


G, 


GT 




.9v 


\9a\ 


(38) 


Gombining these measurements with the bounds on bp in Eq. (30) gives a more precise 


25 
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Figure 2 (Color online) Bounds on left- and right-handed tensor couplings (90% C.L.). The 


measurement of the /3-i^ correlation in ®He (Gliick 


1998 


Johnson et al. 1963) gives a ring-shaped 


boundary. The boundary of measurements of the /3-asymmetry in the pure Gamow-Teller-decay of 
®^Co (Wauters et al. 20101 is given by dashed lines, the measurement only constrains left-handed 
couplings (Eq. (35)). The strongest bounds on left-handed couplings are from measurements of 


the /3-longitudinal polarization Pp/PcT in Eq. (38) (Carnoy et al. 
combined with the constraint on bp. 


1991 Wichers et al. 1987), 


left-handed tensor bound, but it does not constrain right-handed couplings. 


Figure shows the best constraints on tensor couplings. We use the Pf/Pgt values 


(Carnoy et al, 1991 Wichers et ail 1987), the /3-i/ correlation in °He (Gliick, 1998 Johnson 


et al. 1963), and the /3 asymmetry in ®°Co (Wauters et al., 2010) (see Tab. Ill) to hud 
the best bounds for nuclear searches, using minimalization. For the Pf/Pgt values we 


have included the limits on scalar couplings in Eq. (34). The combined £t for real tensor 
couplings gives, at 90% C.L., 


-0.3 X 10“2 < < 0.6 X lO'^ , 

\9a\ 

-6 X 10“2 < < 6 X 10“2 . 

\9a\ 


(39a) 

(39b) 


Reducing the limits will require increased statistics and experimental improvements (Sec. IV.D) 
Further constraints from /3 decay come from mixed decays which we discuss next. 
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3. Tensor constraints from neutron and mirror nuclei 


Mirror transitions are mixed transitions and therefore sensitive to both scalar and tensor 
interactions. Mirror decays might be used to improve the bounds of pure Fermi and Gamow- 
Teller transitions discussed above. At this point only the neutron can be considered. The 
prospects of using mirror nuclei are discussed at the end of this subsection. The neutron can 


serve as a laboratory for studying a range of fundamental interactions (Abele, 2008 Dubbers 


and Schmidt, 2011 Nico, 2009)). In neutron /3 decay, the main focus lies on determining the 
SM parameters and A = Qa/Q v- Non-SM values are included by allowing A to be complex 
and/or by allowing for scalar (A^, Ar) and/or tensor an) interactions. We still consider 
only real couplings, and defer to Sec. V.A.l and Sec. V.A.2 for complex A and scalar and 
tensor couplings, respectively. To clarify the role of possible left- and right-handed scalar 
and tensor contributions, we keep the simplifying assumptions that the V and A couplings 
are those of the SM. For neutron decay, with Mqx = ^/3 and Mp = 1, the ft value is given 
by 


27r^ ln(2) 


1 /ftn — 


1 + 


9 ^ 

9v 




+3A" 1 + 


9t 

9A 


“l + “f? 


-‘’(Da-* 


(40) 


The current value recommended for the lifetime is Tn = 880.3(1.1) s (Olive et al.. 2014), 


which is nearly 6 seconds lower, but with the same error, as the recommended value of 
2008. Of course, this affects the SM values for Vud and A, but cross-checks with other 


correlation coefficients are possible, allowing for consistency of the SM parameters (Wietfeldt 


and Greene, 2011). Including scalar and tensor contributions increases the number of degrees 


of freedom and such cross-checks are no longer possible. The observable ftn is most sensitive 
to ap, because of the partial Gamow-Teller nature of neutron decay. One can combine 
various correlation coefficients from neutron decay to extract A, while allowing for non-SM 
contributions. In combination with the experimental results from the superallowed Fermi 
transitions {bp and iFt), improved bounds on tensor contributions can be obtained. For 


example, with the recent limits on A from UGNA and PERKEOII (Mendenhall et al., 2013 
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Figure 3 (Color online) Contour plot of the 1, 2, and 3 a contours, derived from the selection of 


available data listed in Table III In the fitting procedure we minimized Al, Aji, and A. Notice the 
scale difference of the two axis. 


Mund et al. , 2013) and neglecting right-handed nentrinos {Aji = 0, = 0), it is possible to 


obtain an analytical bonnd on (Pattie et al. , 2013). Allowing for right-handed nentrinos 
reqnires a fitting procednre. 


A complete set of nentron correlation data has been compiled by Dnbbers and Schmidt 


(2011). More recent resnlts are obtained with the PERKEOII setnp (Mnnd et al. 2013) and 


from the UCNA collaboration (Mendenhall et al. , 2013). Combined with the bounds from 


pure Fermi and Gamow-Teller transitions a fit can be made to obtain all relevant parameters 


(A, Ai, Ar, ai, anda^j) in a consistent way. This was recently done by Wauters et al. (2014), 


to extract both left-handed and right-handed tensor-coupling limits. Their fitting method 
entails a grid search. For all ai and values, a value of was obtained by minimizing 
for the other 3 parameters. With this 2D x^ surface a contour plot can be made, by plotting 
the equal = x^ — Xo lines, where Xo is the minimal 

Figure 1^ shows the contour plot for the 1, 2, and 3 a (Ay^ = 1,4, and 9) bounds 
obtained with this method and by using the most relevant experiments listed in Table [TTTj It is 
important to note that the neutron lifetime requires the value of Vud- The most precise value 























for Vud is obtained from the J^t of superallowed decays (Hardy and Towner, 2014), under the 


assumption of no scalar interactions. We have corrected for this by using Eq. (A20) for the 


neutron lifetime. For the neutron lifetime we use the average value of the PDG (Beringer 


et al ., 2012). For the correlation coefficients the averages of the PDG cannot be used, because 


these are obtained by assuming only SM interaction. The possible different dependence on 
the Fierz-interference term is therefore not included. We consider the different values of A 


separately, for which we have calculated the energy dependence with Eq. (26). We have 
included the measurement of H, although for neutron decay this coefficient actually has a 


reduced sensitivity to the Fierz term b and to A, see Eq. (A26). 
We End at 90% G.lfl 


-0.3 X 10-2 < < 0.06 X 10-2 , (41a) 

I^aI 

-4.6 X 10-2 < -^aR < 4.6 x IO -2 , (41b) 

I^aI 

-0.1 X 10-2 < ^Al < 0.3 X 10-2 , (41c) 

9v 

-5 X 10-2 < ^Ar < 6 X 10-2 , (41d) 

9v 

1.2659 < A < 1.2746 . (41e) 


The extracted value of A has a much larger error compared to A = 1.2723(23) from PDG. The 


scalar bounds are the same as the bounds in Eq. (34), but the tensor bounds are improved 
because of the inclusion of the neutron data. Especially the positive bound for aR is reduced 


as compared to Eq. (39). This is caused by the large spread in experimental values for A. 


Using only the two most recent values of the PERKEOII setup (Mund et al.. 2013) and 


from the UGNA collaboration (Mendenhall et al., 2013) gives —0.3 x 10 2 < grpaL/\9A\ < 
0.2 X 10-2. For the tensor bounds, the neutron lifetime has a large influence (Wauters et al 


2014). We therefore anticipate that the error in the neutron lifetime and the spread in A 


will soon give the dominant error on the limit on tensor couplings. 

Recently, also mirror decays have been used to constrain tensor couplings. The strong 
constraint on bp from superallowed Fermi decays, can be combined with measurements on 


mirror nuclei, to derive a value for bar- In Severijns et al. (2008) a complete survey of 

^ Bounds are extracted by scanning the 2D + 1-64^ surface for scalar (Al^r) and tensor (oL.i?), while 
for A we used the ID probability density. 
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Isotope 

Parameter 

Decay 

{nie/Ee) Value 

Error 

^He 

dcT 

/?-, GT 

0.286 

-0.3308 

0.003 

140/iOG 

Pf/Pgt (Eq- (38 

)) (9+ 

0.292 

0.9996 

0.0037 

26m^jy'30 

P Pf/Pgt (Eq. (38 

)) (9+ 

0.216 

1.003 

0.004 

32 Ar 

dp 

/?+,F 

0.191 

0.9989 

0.0065 

38m 

dp 

/3+,F 

0.133 

0.9981 

0.0045 

60Go 

Act 

/?-, GT 

0.704 

-1.027 

0.022 

0+ ^ 0+ 

bp 

/3+,F 

0.2560 

-0.0022 

0.0026 

n 

T (Eq. (A20)) 

/?-, F/GT 

0.655 

880 s. 

0.9 s. 

n 

An 

/?-, F/GT 

0.56 

-0.11952 

0.00110 

n 

An 

/?-, F/GT 

0.534 

-0.11926 

0.00050 

n 

An 

13-, F/GT 

0.582 

-0.1160 

0.0015 

n 

An 

/?-, F/GT 

0.558 

-0.1135 

0.0014 

n 

An 

(3-, F/GT 

0.551 

-0.1146 

0.0019 

n 

Bn 

(3-, F/GT 

0.594 

0.9801 

0.0046 

n 

Bn 

(3-, F/GT 

0.63 

0.9802 

0.0050 

n 

Qjn 

/?-, F/GT 

0.655 

-0.1054 

0.0055 


He and McKellar 


Wichers et al. 


Adelberger et al. 


Wauters et al. 


Hardy and Towner 


Beringer et al. 


Gliick 


(1998) 


(1993) 


Carnoy et al. 


(1991) 


(1987) 


Gorelov et al. 


(1999) 


(2005) 


( 2010 ) 


(2009) 


( 2012 ) 


Mendenhall et al. 


Mund et al. 


Liaud et al. 


(2013) 


(1997) 


(2013) 


Yerozolimsky et al. 


Erozolimskii et al. 


Bopp et al. 


(1997) 


(1991) 


£9g 


Serebrov et al. 


(1998) 


Schumann et al. 


Byrne et al. 


(2007) 


( 2002 ) 


Table III Experimental values used to construct Fig. The values for {nielEe) are mostly not 
calculated by the experimental groups and are derived with Eq. (26), except for the O"*" —)• 0"*“ 


decays, for which we use the value derived in Pattie et al. (2013). Averages from the PDG are only 


used for the r (Beringer et al., 2012), since different measurements of A and B might also have a 


different energy dependence, which is not taken into account in the PDG averages. We have taken 
all experimental values for A used by the PDG. 


J^t values of the available mirror transitions is given. For T = 1/2 transitions the relation 


between the J^t values of the mirror and superallowed 0+ —?■ O’*" is given by (Severijns et al. 


2008) 


2Tt' 


0+-5>0+ 




l + |[.4i + /l|]-27(u) 


0 +^> 0 + 


9v ^ 


1 + I [yli + Al] + [1 7 4^2 + 4„2^] ± 27 (fj) (iJAi - 2£|aip2) 


( 42 ) 
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where Ja/ fv = 1.0143(29) is the ratio of the axial-vector and vector statistical rate functions 
(Severijns et al.\ 2008[). The inverse energy dependence of the superallowed Fermi decays is 


denoted by {me/Ef,)^ and calculated in 
ai can be extracted from 


Pattie et al. 


(2013). If p is known, a value for 


The mirror /?+ decay of ^®Ne to was recently studied to determine the lifetime of 


^®Ne (Broussard et al ., 2014). In this work, the effectiveness of the method described above 
is shown. For mixed decays an independent measurement of p is necessary. For ^®Ne, p = 


1.5995(45) (Calaprice et al.. 1975), was derived from the measurement of the jS asymmetry 


A. Neglecting quadratic couplings in Eq. (42) and using the extracted value Et = 1719.8(13) 


s with (me/Ee) = 0.387022(18) from Broussard et al. (2014) a limit on box is derived. For 


left-handed tensor couplings this gives at 90% C.L. (Broussard et al., 2014) 


1.5 X 10-2 < < 0.12 X 10-2 

19 a\ 


(43) 


The bounds are only an order of magnitude less precise than the combined limits in Eq. (41), 


and show the potential for this kind of measurements for improving the existing bounds. 


4- Tensor constraints from radiative pion /3 decay 


In Bychkov et al. (2009) limits on tensor couplings are derived from radiative pion decay. 


TT^ —)■ e"' 


■h'e + 'y. These bounds can be translated into bounds on ax (Bhattacharya et al. 


2012) by using estimates for the pion form factor (Mateu and Portoles, 2007). Assuming no 


right-handed couplings and using gx = 1.047(61), a limit at 90% C.L. is found. 


- 1.9 X 10-^ < < 2.3 x 10-^ . (44) 

\9a\ 

These bounds are the strongest bounds on tensor couplings from a single decay experiment 
and show that future /3-decay experiments should probe ax < 10-^ and beyond, in order to 
improve these existing limits. 
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5. Pseudoscalar constraints 


Pseudoscalar interactions have so far been neglected in /3-decay searches, since they are 
strongly suppressed because the nuclei are nonrelativistic. The suppression of these terms is 
0{1/M), where M is the nucleon mass. However, in (3 decay, the pseudoscalar interactions 


are always multiplied by gp, the pseudoscalar form factor discussed in Eq. (22). The large 


value gp = 349(9) (Gonzalez-Alonso and Camalich, 2014) largely cancels this suppression, 
and /3-decay experiments might be used to probe these interactions. There are, however. 


already strong constraints on pseudoscalar couplings from pion decay (Bhattacharya et al. 
[Herczeil [200l| ). 


The ratio = r(7r —)■ eh')/T{7i /ii/) is sensitive to pseudoscalar couplings dehned by 


C = 


GpVud 


e(l-75)i/e + A^e(l+75)z/e M75d , 


P 7^.1 


(45) 


where we have neglected flavor-changing couplings, which can be found in Bhattacharya et al. 


(2012). The ratio where R^^ is the measured value, is sensitive to electron and 

muon pseudoscalar couplings, and respectively. If these couplings are such that 

A^^^^/nie = A^G)I their contributions to the ratio cancel and no bounds on pseudoscalar 
interactions can be obtained. Since there is no reason to assume such a cancellation, we can 
place bounds on pseudoscalar interactions, because these would show up as Rt,/R^ ^ 1. 


The current best value for this ratio is R.,^/R^ = 0.996(3) (Beringer et al, 2012 Cirigliano 


and Rosell, 2007), which leads to (90% C.L.) (Bhattacharya et al., 2012; Cirigliano et al. 


2013b) 


-1.4 X 10“^ < < 5.5 X 10“^ , (46a) 

-2.8 X 10“^ < 71^ < 2.8 X 10“^ . (46b) 


In /3 decay the pseudoscalar term shows up in Gamow-Teller and mixed decays. The 
most relevant to experiments are its contributions to the Fierz interference term. 


9t _ , r, fi'P ^pEo — Ef, 


bcT = ± 2y—^A^ 


'\9a\ 


\9a[ 


M 


(47) 


which enters with the usual {me/Ef.) suppression. The {Eq — Eg)/M term is responsible 
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for the suppression of pseudoscalar contributions, however, because gp{Eo — Ee)/M ~ 0.4 
pseudoscalar interactions are still suppressed compared to tensor interactions. Given the 


current limit on ap, improving the bounds in Eq. (46a) seems unlikely in the near future. 


The pseudoscalar couplings in Eq. (46) can also be translated into bounds on scalar and 
tensor couplings. If scalar and tensor interactions are present at the new physics scale A, 
they will mix via radiative loop corrections, and pseudoscalar couplings will radiatively be 


generated (Campbell and Maybury, 2005 Herczeg, 1994). Current limits are at the level of 


Bhattacharya et al. ^ 

2012 

Cirigliano et al. 

I2013alb 


\Al\ < 8 X 10“2 and |Ah| < 5 x 10“^ , (48a) 

la^l ^ 2 X 10“^ and |a_R| < 1.2 x 10“^ , (48b) 


and depend logarithmically on the scale of new physics A, for which A = 10 TeV is used. 


These bounds are of the same order of magnitude as global-ht limits from jS decay in Eq. (41), 
except for the bound on ap, which is an order of magnitude better. However, because the 


constraints for right-handed currents rely on the flavor structure of new physics (Cirigliano 


et al. ^ 2013b), we do not further consider these bounds. 


6. Left-handed scalar versus tensor 


In Sec. IV.C we discuss exotic couplings involving right-handed neutrinos. If right-handed 
neutrinos are absent, or too heavy to be energetically allowed in jS decay, right-handed 
neutrino couplings, i.e. Ar and an, can be neglected. The resulting reduction of parameter 
space allows us to use mixed decays to £t the correlations between left-handed tensor and 
scalar couplings. Figure shows these correlations. For the complete set of data listed in 
Table [Till we hnd at 90% C.L. 


-0.1 X 10-2 < ^Al < 0.3 X 10-2 , 

9v 

(49a) 

0.2 X 10-2 < < 0.06 x IO -2 , 

\9a\ 

(49b) 

1.2715 < |A| < 1.2744 . 

(49c) 
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Figure 4 (Color online) Contour plot of the 1, 2, and 3 a contours, derived from the selection of 


available data listed in Table III without right-handed couplings, i.e. = or = 0. 


These bounds are not significantly different from the bounds from the complete £t in 


Eq. (41). For comparison: limits on right-handed couplings from neutron decay alone are 


found in Konrad et al. (2010) and Dubbers and Schmidt (2011). 


B. Constraints from LHC experiments 


Low-energy experiments are mostly viewed as complementary to high-energy collider 
searches for BSM physics. Experiments at LHC can place bounds on new physics by looking 


for the on-shell production of new particles, as done in searches for a Wr boson (Eq. (13)) 


or supersymmetric particles. We focus here on the effect of a Wr boson, because this has 


been studied complementary by precision decay experiments and by LHC, e.g. Dekens and 


Boer (2014b). At the LHC, Wr is searched for by considering its possible decay channels. 


In the Wr —?■ tb channel, such direct searches at the CMS experiment constrain Mr > 2 TeV 


(Chatrchyan et al, 2014). Constraints from the Wr —)■ eu channel are similar, but depend 


on assumptions for the right-handed neutrino. Constraints from neutral-kaon mixing give 


Mr > 3 TeV (Bertolini et al, 2014). 


In (3 decay, strong limits come from CKM unitarity tests, for which the best bound is 
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(Hardy and Towner, 2014| ) 


\Vud\" + iK^r + iKfti" = 1.00008(56) , 


(50) 


which uses the value of 14® from 

Moulson 

c 

2013 

). 

The error has equal contributions from 

Vud and Vus- Following 

Hardy and Towner 

(2009 

), this leads to a constraint on rrr, i.e. 


left-handed lepton couplings and right-handed quark couplings, of 

- 4 X 10“^ <aLR<5x 10“^ , (51) 


at 90% C.L. The precision of both Vud and Vus should improve simultaneously for such a 
test to remain signihcant. 


In (3 decay, some correlation coefficients are sensitive to aLR-, aRL, and aRR, where the lat¬ 
ter two are only present if light right-handed neutrinos are assumed. For example, the mea¬ 


surements of Pr/Pgt (Carnoy et al., 1991, Wichers et al. 


and Act in ®°Co (iWauters 


et al., 2010), are used to constrain parameters of manifest LR-symmetric models. Such 
models have a P symmetry, such that for the CKM matrices = ±1^^. There is no 
additional spontaneous CP violation, so a; = 0. In this simplihed model, rrr = PaRR ~ ^ 

and rrr = 6 = {M 1 /M 2 Y. Measurements of Pf/Pgt limit the combination 6 ■ ^ and do 
not give additional bounds, because of the strong bound on ^ from unitarity tests given in 


Eq. (51). Because ^ is strongly constrained, /3-decay experiments can only constrain rrr and 


thus the mass of the Wr. Derived limits are of the order of 200 GeV (Gorelov et al ., 2005 


Wauters et al.. 2010), an order of magnitude below the bound from the LHG experiments 


presented above. In fact, when assuming manifest LR symmetry, the strongest bound on 


Wr comes from the Kr-Ks mass difference, from which Wr > 20 TeV was derived (Maiezza 


and Nemevsek, 2014). 


Besides constraining new physics by searching for direct on-shell production, it is also 
possible for colliders to constrain exotic couplings. When the mass of the non-SM particle 
exceeds the energy accessible at LHG, the new particles cannot be produced on-shell, but 
their effects can still be found in deviations from the SM predictions. In that way, the exotic 
interactions in Eq. will also manifest themselves in proton-proton collisions. This makes 
it possible for LHG data to constrain the same tensor and scalar couplings relevant in (3 decay 
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(Bhattacharya et al., 2012 Cirigliano et al ., j2013b ). 

In particular, the pp —)■ e + MET + X channel is considered, where MET signihes Missing 
Transverse Energy. This channel is closely related to (3 decay, since it involves the ud —)■ eh 
process at quark level. At the LHC, both the ATLAS and CMS detectors are used to search 


for new physics in this channel (Aad et al ., 2012 Chatrchyan et al., 2012), by searching for 


an excess of events predicted at a large lepton transverse mass cut friT. At large friT, the SM 
cross section approaches zero more rapidly than the cross sections for new physics, making 
the sensitivity to non-SM physics larger at high momenta. The total cross section is 


{rriT > ttit) — o'sm(1 + 

+as{\AL\^ + \An\^) + lM\aL\^ + M^^ 


(52) 


where usm is the SM cross section and (Jrst are the cross sections for new physics. The ex¬ 


plicit form of (JsM and is given, to lowest order in QCD corrections, in Cirigliano et al 


(2013b). The coefficients rlr and url cannot be constrained, because their contribution is 


proportional to ctsm, and therefore small at large rriT. 

With the expected number of background events and the number of actual observed 


events, one can place an upper limit on the number of new physics events, (Bhattacharya 


et al, 2012). This translates into an upper limit for a, and finally into bounds on exotic 


couplings. First bounds were derived by Bhattacharya et al. (|2012|), updated bounds are 
given in 


Naviliat-Cuncic and Gonzalez-Alonso (2013) 


The bounds are derived by using the experimental data of Khachatryan et al. (2014) 
at an integrated luminosity of 20 fb“^ and at a center-of-mass energy of = 8 TeV. 


Naviliat-Cuncic and Gonzalez-Alonso (2013) also gives the combined limits for scalar and 


tensor couplings, assuming only left-handed couplings. In Table W we give the 90% G.L. 
bounds, obtained by allowing one exotic interaction and putting all other couplings to zero. 
To compare these results with /3-decay constraints, we use the values from the global fit 


in Eq. (41) and the form factors gs = 1.02(11) (Gonzalez-Alonso and Gamalich, 2014) and 


gx = 1.047(61) (Bhattacharya et al., 2014). Because the errors on the form factors are not 


Gaussian, we use the R-fit method described in Bhattacharya et al. (2012), which treats all 


the values in a 1 a interval with equal probability. Therefore, only the lower bounds are 
important. We stress again that the reduction of the error in gs and gx is important to 
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I^lI 

\^r\ 

lei'll 


/3 decay 

2.5 X 10-3 

6 X 10-2 

3 X 10-3 

4.6 X 10-2 

LHC 

6 X 10-3 

6 X 10-3 

2 X 10-3 

2 X 10-3 

Neutrino 

- 

1 X 10-3 

- 

1 X 10-3 


Table IV Comparison between /3-decay limits on left- and right-handed scalar Al and and 


tensor couplings ai and ur, constraints from LHC data (Naviliat-Cuncic and Gonzalez-Alonso 


20131, and from the neutrino mass (Ito and Prezeau 2005|. Constraints are at 90% C.L., and all 


couplings are assumed to be real. 


make meaningful comparisons between the different experiments. 


Table IV shows that the LHC constraints on left-handed couplings are comparable to 
/3-decay constraints, while for right-handed couplings the LHC constraints are an order of 
magnitude better than the /3-decay limits. The current status is illustrated in Fig. and 
Fig. [Navihat-Cuncic and Gonzalez-Alonso (2013) also make a projection for the 14 TeV 
run at 50 fb“^, and hnd that the expected bounds are a factor 3 better. 


C. Neutrino-mass implications 


Besides strong bounds from LHC experiments on right-handed interactions, there are 
also bounds from the neutrino mass. In the SM, neutrinos are assumed to be massless, but 
neutrino oscillations indicate the existence of at least two massive neutrinos. A direct upper 
limit on the neutrino mass comes from the shift of the end-point of the /3 spectrum. Recent 


measurements of the /3 spectrum of ^H give rriy < 2 eV (95% C.L.) (Aseev et al, 2011 


Kraus et al, 2005). The KATRIN experiment aims to improve these limits by an order of 


magnitude (Otten and Weinheimer, 2008|). Other bounds on the neutrino mass are derived 


from cosmological observations; WMAP (Hinshaw et al, 2013) limits J^mi, < 0.44 eV and 


a recent study of Planck (Ade et al, 2014), in which Planck data is combined with neutrino 


oscillation data, gives a similar limit < 0.15 eV, for three degenerate neutrinos. 

In Eq. the couplings aRR, aRi, Ar, and aR involve right-handed neutrinos. These 
couplings can only be generated if the decay to right-handed neutrinos is kinematically 
allowed, i.e. if right-handed neutrinos are light enough to be created in the decay. The 
possibility of these light right-handed neutrinos has been considered in various new physics 
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(a) Constrains and an 


(b) Constrains qrl 


Figure 5 The two-loop contribution to the neutrino mass, where the box indicates the exotic 


couplings. The cross indicates a mass insertion, with (a) mq = 4 MeV (b) rUe = 0.511 MeV (Ito 


and Prezeau 20051. For Majorana neutrinos one can substitute ur —)■ u' 


L- 


scenarios as a possible dark-matter candidate. If right-handed nentrinos are very heavy, as 
is snggested in many see-saw mechanisms, we can omit all exotic conplings with hrst index 

R. 


Prezean and Knrylov (2005) showed that the small nentrino mass also limits the presence 


of exotic conplings in low-energy experiments that involve a (light) right-handed nentrino 


(Prezean and Knrylov, 2005). For /? decay this strongly constrains the conplings and 


o-RL (Ito and Prezean, 2005). Nentrino masses can be either Dirac or Majorana 


where z/£ = * 7270 ^'^, or a combination of the two. However, the following resnlts 
are general and apply to both types. Conplings to right-handed nentrinos contribute to the 
neutrino mass via loop interactions. Figure]^ shows the leading two-loop contribution to the 
neutrino mass, where the box indicates the non-SM coupling to right-handed particles. The 
cross indicates the mass insertion needed to couple two fermions with different chiralities. 
Here, the chirality-changing interactions are either proportional to (a) the quark or (b) 
the electron mass. In a power-counting scheme, one-loop contributions are in general less 
suppressed than two-loop contributions. However, the two-loop diagrams in Fig. are 
enhanced by the IF-boson mass, while the one-loop diagrams are only suppressed by the 
light-fermion mass. This makes the two-loop contribution dominant, as the additional loop 
suppression of l/(47r)^ is diminished by the heavy IK-boson mass. 


One can estimate the two-loop contribution to the neutrino mass by considering only the 
logarithmic part of Fig. The analytic parts are renormalization-scheme dependent and 


are therefore neglected (Prezeau and Knrylov, 2005). By using dimensional regularization 
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the contribution to dmy is estimated as (|Ito and Prezeau 2005) 




2 ^ fi- 


(dvr) 


In 




( 53 ) 


where a = «_r, orl} are the exotic couplings from Eq. 0 , g = 0.64 is the 

gauge coupling, ruf is the inserted fermion mass, and fi is the renormalization scale, which 
should exceed the heaviest mass in the interaction, /i > rrit, where m* is the top-quark mass. 
Assuming that the loop corrections do not exceed the mass of the neutrinc]^ i.e. drrii, < 


setting rrig = 4 MeV, g, = I TeV, and rrit, < 0.15 eV in Eq. (53) gives 


I^rlI ^ 10 ^ , 

|Arr| , |Arr| , I^rI ^ 10 ^. 


(54a) 

(54b) 


In Table IV we compare these limits with current right-handed /3-decay bounds and bounds 
from LHC. The estimates from the neutrino mass are currently the strongest bounds on 
right-handed currents. They are more than an order of magnitude stronger than the /3-decay 


bounds, and comparable to the LHC bounds. For the bounds in Eq. (54) we have used the 


updated neutrino mass from the Planck space observatory, which might further improve in 
the future. The given bounds are conservative estimates, but nevertheless they show the 
large impact of the neutrino mass on /3-decay measurements. Even stronger constraints of 
0(10“^) from the neutrino mass have been derived in the unpublished thesis of Wang (2007). 


D. Conclusions and outlook 

We summarized the current status of the bounds on real right-handed vector, scalar, 
pseudoscalar, and tensor interactions in /3 decay. We compared these bounds with those 
obtained from proton-proton collisions at the LHC experiments and the upper limit on the 
neutrino mass, mainly focusing on scalar and tensor interaction. The best current bounds 


are given in Table IV We distinguished between bounds on left- and right-handed scalar and 


tensor interactions, where left or right denotes the chirality of the neutrino. The constraints 
on left-handed interactions are equally constrained by the LHC and /3-decay experiments. 


^ There might be scenarios in which this is not obeyed, but these scenarios would have to be fine-tuned. 
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8v 

Figure 6 (Color online) Scalar bounds from nuclear /3 decay as in Fig. [^combined with limits de¬ 
rived from the neutrino mass (horizontal line) and constraints from the LHC experiments (circular 
bound). 
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Figure 7 (Color online) Tensor bounds from nuclear /3 decay as in Fig. [^combined with limits de¬ 
rived from the neutrino mass (horizontal line) and constraints from the LHC experiments (circular 
bound). 


On the other hand, /5-decay experiments measuring right-handed interactions would have 
to improve orders of magnitude to compete with the bounds from the LHC experiments 
and the neutrino mass. This is illustrated in Fig. for scalar interactions and in Fig. 
for tensor interactions. Table |V] projects the competitive accuracies required for different 
/5-decay parameters. For left-handed currents we give the necessary precision to compete 
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with projected future LHC bounds (Naviliat-Cuncic and Gonzalez-Alonso, 2013). For right- 
handed bounds, we give two accuracies. The hrst corresponds to the required sensitivity 
to compete with current LHC bounds, the number in brackets corresponds to the required 


precision to compete with the bounds from the neutrino mass (see Table IV). 

The bounds on left-handed couplings are best pursued via measurements of the Fierz 
interference coefficient b. For left-handed scalar couplings Al the bound is most stringent 
because of the vast effort in the study of super-allowed Fermi transitions. These studies also 
provide the best current value for Vud- The left-handed tensor coupling ai requires a larger 
effort, for which several measurements need to be combined. The best current bounds are 
from the global £t in which neutron and nuclear data are combined. In this fit, especially the 
uncertainties in the neutron lifetime and the A coefficient of the neutron have a signihcant 
impact. We pointed out that the large spread in the available A measurements influences 
the obtained bound signihcantly. The Gamow-Teller part Bgt of the Fierz interference term 
and Vud can also be constrained in mirror nuclei, in analogy to the superallowed Fermi tran¬ 
sitions. However, this also requires the measurement of at least one correlation coefficient. 


Measurements with this aim are undertaken (Ban et al.. 2013). 

In Gamow-Teller transitions, measurements of the Fierz interference term Bgt allow for 
bounds on the left-handed tensor terms. In Seattle, a ®He factory has been set up to study 


this term. The lifetime of °He was already measured with high precision (Knecht et al. 


2012 ), but the shell-model calculations are not sufficiently accurate as yet to search for 


tensor interactions. One straightforward, but not so simple, approach is to measure the 
decay spectrum precisely. This would give access to box- These measurements would also 
have to consider contributions from the SM weak-magnetism (cf. Eq. 0). Measurements 
of bcT from electron-antineutrino correlation dpu and the spectrum are both ongoing and 


being set up (Aviv et al., 2012 Flechard et al, 2011, 2008 Knecht et al, 2011 Naviliat- 


Guncic, 2014 Severijns, 2014). If these measurements reach 6 < 10 they would allow for 


a strong limit on a^. Such a precision is necessary to compete with the projected bounds 
from the 14 TeV run of the LHG. In neutron decay, many efforts are undertaken to improve 


the measurements of a^i, and A (Baessler et al, 2008, 2014, Konrad et al, 2012 Markisch 


et al.. 2009 Pocanic et al ., j2009 Wietfeldt et al, 2009). For comparison, limits on the Fierz 


terms from neutron decay alone are found in Konrad et al. (2010) and Dubbers and Schmidt 


(2011), including limits derived from the electron energy dependence of the /^-asymmetry 
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Axp(-E) alone. 


Right-handed interactions, which imply the existence of a light right-handed nentrino, 
do not interfere with the SM interactions and can therefore only be measured directly, i.e. 
via quadratic terms. This makes it difficult to reach the sensitivity obtained for left-handed 
couplings. In 13 decay, the right-handed tensor coupling can be constrained by measuring 
the /3-z/ correlation, djs^. The best measurement in pure Gamow-Teller decays of stems 


from the measurement in “He (Johnson et al., 1963). Many efforts are undertaken to improve 


this limit in ®He (Aviv et al. ^ 2012 Couratin et al., 2012 Knecht et al., 2011). A dedicated 
effort to limit right-handed tensor couplings is ongoing in ®Li, for which the daughter nucleus 
®Be breaks up into two a particles, ®Li —)■ e~ + d + 2a. The acr coefficient can be measured 
by measuring the (3-a correlation, and by taking advantage of the increased sensitivity due 
to the population of a state in ®Be. After putting the Fierz term 6 = 0, such that only 
right-handed interactions are constrained (|Li et al., 2013), one finds 


^|aR|<8xlO-^ (55) 

\9a\ 

The bound reaches the precision of the combined fits, but when considering the LHC or neu¬ 
trino bounds the experiment would have to improve by more than three orders of magnitude 
to compete (see Table 0, 

When comparing tensor and scalar bounds from different fields, the form factors gs and 
qt are important. Lattice QCD calculations have made enormous progress, and will continue 
to do so in the next period. The lattice prediction of qa will hopefully reach the experimental 
precision soon, which would allow for a cross-check between the experimental value and the 
theoretical lattice value. 

Besides scalar and tensor searches, we also discussed searches for H -|- A and pseudoscalar 
interactions. Pseudoscalar interactions are less suppressed than previously thought, due to 
the large value of gp. However, strong bounds exists from radiative pion decay, and pseu¬ 
doscalar interactions can still be neglected in the upcoming /3-decay experiments. Strong 
constraints on V + A currents are extracted from CKM unitarity tests, to which /3-decay 
experiments contribute by providing the most accurate value of Vud- Besides this, measure¬ 
ments of correlation coefficients can be used to constrain parameters of (manifest) left-right 
symmetric models. For these specific models, strong limits from the LHC experiments and 
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Parameter 

Bound 

Constraint at 90% C.L. 

hcT 

10-3 

ttL < 3 X lO-'^ 

bp 

10-3 

Ai < 5 X 10-^ 

acT 

10-^ (5 X 10-®) 

aR<Qx 10-3 [ur < 10-3) 

dp 

8 X 10-® (2 X 10-®) 

< 2 X 10-3 {Ar < 10-3) 


Table V Required experimental precision on /S-decay parameters to remain competitive with LHC 
bounds, cf. Naviliat-Cuncic and Gonzalez-Alonso (20131. Only the Fermi (F) and Gamow-Teller 
(GT) parts of the Fierz-interference term b and the /3-i/-correlation a are listed. The third column 
gives the corresponding limit on scalar couplings and Ar and tensor couplings ol and aR. 
The Fierz term is the leading term in most /3-correlation experiments (Sec. IV.A). The indicated 
bounds for b assumes that future LHC data lead to bounds indicated in the last column. The a 
parameter is the most direct way to obtain a bound on right-handed couplings, which should be 
the motivation to measure a. Here the current bounds of thr LHC are assumed, while the value in 
parentheses is the required accuracy when the bound derived from the limit of the neutrino mass 
is considered (Table [T^. 


the neutral-kaon mass difference exist. Therefore, the signihcance of (3 experiments in these 
experiments is limited to specihc models. 


V. LIMITS ON TIME-REVERSAL VIOLATION 


So far we have only considered the real parts of the exotic couplings. In this Section we 
focus on their imaginary parts. A nonzero measurement of an imaginary coupling would 
imply that time-reversal (T) symmetry and, by the CPT-theorem, CP symmetry is violatecQ 
Becaues of the matter-antimatter asymmetry of the universe, new sources of CP violation 


are expected (Sakharov, 1967). Many models of BSM physics predict such additional sources 


of CP violation, see e.g. Branco et al. (2012); Dekens and Boer (2014a); and Ibrahim and 


Nath (2008). This makes T or CP violation one of the main portals to search for new 


physics. These searches range from experiments at the LHC to atomic-physics experiments. 
As such the observables can be quite diverse. With advances in theory, in particular via EFT 


methods, relations between the different observables have become more clear (cf. Sec. IV.B 


and Sec. IV.C). 


In any Lorentz-symmetric local field theory, CP violation is equivalent to T violation, according to the 


CPT theorem. For CPT-violation, see Sec. VI 
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In this section we focus on the connection between T-violating observables in f3 decay and 
the bounds on electric dipole moments (EDMs). The P- and T-odd EDM measurements 


are a powerful probe of CP violation beyond the SM (Pospelov and Ritz, 2005). High- 
precision EDM searches have been made for the neutron, paramagnetic and diamagnetic 
atoms, and molecules. The EDM is a static observable, and, therefore, allows for very precise 
atomic-physics experiments. It is also a background-free observable, because the electroweak 
SM contributions to the EDM are strongly suppressed. Therefore, EDM experiments give 
strong limits on new T-violating physics. BSM physics contributions to the EDM can be 
parametrized by dimension-6 operators (de Vries et al ., 2011a[b i 2013, 2011c). At low energy 
this leads to a relation between the T-violating correlations in (3 decay and EDMs. 

Many correlation coefficients in (3 decay depend on the square of the underlying coupling 
constants. As such they depend only on the imaginary couplings squared, which are therefore 
difficult to access. A more direct way to probe imaginary couplings is to consider the T-odd 


triple correlations J-{pe^Pv) and (Tg-(</xpe) multiplied by the D (Eq. (17)) and i? (Eq. (18)) 
coefficients, respectively. The hrst is P-even and T-odd, while the latter is P- and T-odd. 
They probe left-handed imaginary couplings, which are absent in the SM. 

Since the interactions contributing to D, R, and EDMs are generated by the same oper¬ 
ators, a limit on the EDM also limits the D and R coefficients. We consider these relations 
and discuss the relative precision of the two types of experiments. 


A. Limits on triple-correlation coefficients in f3 decay 

A hnite D coefficient arises from the interference between the imaginary parts of the 
left-handed vector couplings and is proportional to Im The R coefficient arises from 
the interference between the imaginary parts of scalar or tensor couplings and SM couplings, 
making this coefficient sensitive to both Im and Im an. 

The SM contributes to both the R and D coefficients through electromagnetic hnal-state- 
interactions (ESI) and through SM CP violation. The ESI are only motion-reversal odd, 
i.e. the initial and hnal state are no longer interchangeable, due to radiative corrections. In 
this way, ESI mimic time-reversal violation, but in fact are T-even. We will denote their 
contributions by R/ and Df, and write D = Dt + Dj and R = Rt + Rf (Herczeg, 2005), 
where Dt and Rt are the true T-violating contributions. The contributions from ESI are 


(Herczeg, 2005 
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comparable to the current experimental precision and depend on the momentum of the f3 
particle. We will discuss their values for specihc isotopes later. True T violation in the SM 
arises from the CP-violating phase of the CKM matrix and the QCD 6*-term. These sources 


only contribute at the level of (9(10 (Herczeg and Khriplovich, 1997), much below the 
current experimental precision. 


1. D coefficient 


To hrst order in exotic couplings, the Dt coefficient can be expressed as (Jackson et al. 

Dt = aDlm aiR , (56) 


1957a) 


from Eq. (A15), with 


a_D — 




(57) 


l + p2 

The D coefficient can only be accessed in mixed transitions, and has been measured in both 
neutron and ^®Ne decay, which have an = 0.87 and ajy = —1.03, respectively. For ^®Ne 


the best measurement is D = 1(6) x 10 ^ (Hallin et al., 1984), and from neutron decay 
D = —0.94(2.10) X 10“^ (Chupp et al.. 2012 Mumm et al.. 2011). 

The value of the ESI depends on the kinematics of the experiment. For ^®Ne the ESI have 


been derived by Callan and Treiman (1967) as Dj = 2.6 x 10 which is of the same 

order as the experimental precision. For neutron decay the ESI were also calculated in chiral 


perturbation theory by Ando et al. (2009). Their derivation reproduces the original result of 


Callan and Treiman (1967). However, Ando et al. (2009) include higher-order corrections. 


which are of order (9(10 ^), allowing for an accurate expression for the ESI, 


^rn.3,x 

Df = ( 0 . 228 ^^ + 1 . 083 ^ 


X 10 


^max 

- 5.88^^ X 10" 


(58) 


Pe Pe Pe 

where the first two terms are the Callan and Treiman (1967[) terms, and the last term 


represents the higher-order corrections. Equation (58) is accurate to better than 1%. For 
the current best neutron experiment the ESI are estimated at i9/ ~ 1.2 x 10-^ dChuppI 


et al, 2012). The uncertainty in ZJj stems from the uncertainty of the (3 momentum in the 
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experiment. The T-violating part of the neutron D measurement gives at 90% C.L. 

I A| < 4 X 10“^ , (59) 

and with qd = 0.87, 

|Im ulrI < 4 X 10“"^ . (60) 

Given the current experimental precision, it is clear that the FSI become increasingly 
more important. In this respect, neutron experiments are favored over nuclei, because the 
FSI can be calculated with a higher precision. Eventually the accuracy to which the FSI are 
known will limit measurements of true T violation. 


2. R coefficient 


Neglecting quadratic non-SM couplings, the Rt coefficient is given by (Jackson et al. 

{an =F ^d) 


1957a) 


Rt = 


an 


-gxlm. aL - - —S'sim , 

\9a\ 2gv 


(61) 


from Eq. (A16), where the upper (lower) sign is for ff (/!+) decay, an is given in Eq. (57), 
and 

bo = , (62) 

with Xjij as given in Appendix]^ The R coefficient can be measured in both mixed or pure 
Gamow-Teller transitions, where the latter limits Im a^. The leading contributions to the 


FSI are given by the Goulomb corrections calculated by Jackson et al. (1957b), 


Zarue, , , . 
Rf — -(T^d + od) ■ 

^Pe 


(63) 


The R coefficient has been measured in the pure Gamow-Teller decay of ®Li, where an = 0 


and bo = 4/3. The FSI give 7?/ ~ 7 x 10 leading to Rt = (0.9 ±2.2) x 10 ^ (Huber et al. 


2003). This constrains at 90% G.L. 


^f^llm < 3 X 10 


-3 


(64) 
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The best measurement of i? in a mixed decay has been obtained for neutron decay, for which 
an = 0.87 and bn = 2.2. Kozela et al.\ (2012) hnd R = (4 ± 12 ± 5) x 10“^. The FSI are 


calculated with Eq. (63). By using the energy distribution seen by the experimental setup 
one obtains i?/ ~ 6 x 10“^ (Kozela et al. ^ 2012). The error in Rf is less than 10%. Rf can 


be neglected given the current experimental precision. At 90% C.L. 


— l.lgf'rini ai — 0.445f5lm Al < 2.4 x 10 


-2 


(65) 


With the constraint given in Eq. (64) one hnds at 90% C.L 


^(^llm Al\ < 6 X 10 


-2 


( 66 ) 


3. Alternative eorrelations 

The measurement of the D coefficient requires the detection of the recoiling nucleus 
instead of detecting the neutrino. This imposes strong experimental constraints on any 
measurement scheme. Current schemes consider atomic trapping in a magneto-optical trap, 
which has led to the best value for the fi-v correlation a. Measuring D requires a modihcation 
of this trap technique, to allow for a polarized sample. It will be extremely challenging 
to achieve high statistical precision and systematical accuracy with this technique. An 


alternative lies in the /3-7 correlations of polarized nuclei (Curtis and Lewis, 1957 Morita 


and Morita 


1957), where the photon with momentum k is emitted from the state populated 


by the (3 decay. In this way one measures the correlation proportional to 


EJ ■ {pe X k){J ■ k) , 


(67) 


when the emission is due to an El transition. The correlation coefficient E oc Im qlr 


is nonzero only for mixed decays. Young et al. (1995) have identihed as a promising 


candidate for such a measurement, since this isotope allows for the comparison between a 
mixed and a Gamow-Teller transition. The latter is insensitive to T violation and can be 
used to test the experimental setup and reduce systematic errors. Secondary beams of high 
intensity can be produced, stopped, and polarized in a buffer gas allowing to measure /d-y 


correlations (Muller et al., 2013) with high precision. Correlations alternative to measuring 
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EDM 

e cm (90% C.L.) 

Reference 

Gonnection to (3 decay 

n 

2.9 X 10-26 

Baker et al. 


2006) 


D 

199Hg 

2.6 X 10-29 

Griffith et al. 

(2001 

)) 

D, R 

205qq 

0.9 X 10-24 

Regan et al. 

( 

2002 


R 

YbF 

141 < 10.5 X 10-28 

Hudson et al. 

(201 

1 ) 

R 

ThO 

\de\ < 8.7 X 10-29 

Baron et al. 

( 

2014) 


R 


Table VI The current best EDM limits of the neutron, diamagnetic Hg, paramagnetic Tl, and 
molecular YbF and ThO. The neutron EDM and Hg can be connected to the D coefficient (and E 
coefficient). Other EDM measurements, except the neutron, can be connected to the R coefficient. 
The limit from molecular YbF and ThO are expressed as a constraint on the electron EDM de- 


R are also possible (the L and M coefficients (Ebel and Feldman, 1957 Jackson et al 


1957a)) but, similar to R, will always require to measure the polarization of the fd particle, 
which is an inefficient process. 

In radiative f3 decay, it is possible to have triple-correlation coefficients without nuclear 
or electron spin (Braguta et al.\ 2002 Gardner and He, 2012, 2013), such as 


Kk ■ (p„ X Pc) . 


( 68 ) 


This coefficient has not been measured, but Dekens and Vos (2015) showed that EDMs 


provide extremely strong constraints on the coefficient K. 


B. EDM limits 


Limits exist for the neutron EDM, the electron EDM, and several atomic EDMs. The 


best current bounds are listed in Table VI where the limits from molecular YbF and ThO 


are expressed as a limit on the electron EDM de- The last column of Table IVT] indicates if a 


connection to the triple-correlation coefficients D and R exists (Khriplovich, 1991 Ng and 


Tulin, 2012). 
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1. Limits on D from EDM limits 


Any new vector interaction that contributes to Ini aLR (and thus to Dt) also contributes 


to nuclear EDMs (Herczeg, 2005 Ng and Tulin, 2012). This makes it possible to translate 


bounds on the EDMs of the neutron and diamagnetic atoms into bounds on Im ulr- The D 
coefficient is P-even and T-odd, while the EDM is both P- and T-odd. Nevertheless, loop 
corrections, containing the W boson, allow for a relation between these observables. 


The relevant CP-odd dimension-six operator is (Ng and Tulin, 2012) 


= j^nR-(“dR 'f'iDRip + h.c. 


(69) 


where c is a complex coefficient, A is the scale of new physics, is the covariant derivative, 
and (p is the Higgs doublet with where is the antisymmetric tensor. Fig. 

shows the energy evolution of this operator. First, electroweak symmetry breaking generates 
the coupling of the W boson to right-handed quarks, 

URrdRW; + C-dRYURW;) , (70) 

where cp acquired its vacuum expectation value v/\/2 and g is the SU{2)l coupling constant. 
The W boson can couple to a lepton current or a quark current. At lower energy, the W 
boson is integrated out. This generates a P- and T-odd four-quark coupling and the lepton- 
quark coupling qlr in (3 decay. The effective Lagrangian is 


= -h {liRl^dReLl^VcR + V„dUR-y'‘dRdn^ui'j + h.c. . 


( 71 ) 


which shows that the two couplings c and rlr have a common origin. They are related by 


Im Rlr = 


Im c 


2V2GfA^ ■ 


(72) 


When evolving to the QCD scale, the second term in Eq. ( [7l| is affected by QCD renormal¬ 
ization. However, this only has a small numerical effect ( [Dekens and Boct 2014a), which 
can be neglected given the uncertainties coming from the calculation of the neutron EDM. 

Bounds on Im c thus lead to an upper limit on Im The dependence of the EDM on 
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Figure 8 Generation of the four-fermion operators that contribute to the EDM (left) and /3-decay 
(right). The boxes denote the four-fermion couplings, c and aiR respectively. The coupling of the 
W boson to the right-handed quarks is generated by the dimension-6 operator in Eq. (691. 


Im c involves theoretical calculations at different energy scales. Especially for diamagnetic 
atoms such as differences in nuclear calculations lead to a large uncertainty in the 

interpretation of the bounds on atomic EDMs. Therefore, we do not consider bounds from 


199 


Hg. No such problem occurs for the neutron, and Seng et al. (2014) and de Vries et al. 


(2013) estimated the link between the neutron EDM and Im c as 


= — 1 X 10 


-20 


Im c 


2V2GfA^ 


e cm . 


(73) 


This result differs by an order of magnitude from the result used in Ng and Tulin (2012), 


which was obtained from An et al. (2010) and He and McKellar (1993). In Seng et al. (2014) 


and de Vries et al. (2013) it was pointed out that, due to use of a relativistic meson-nucleon 


held theory. An et al. (2010) and He and McKellar (1993) overestimate the neutron EDM 


by an order of magnitude. 


The current bound on the neutron EDM \dn\ < 2.9 x 10 e cm (Baker et al ., 2006) and 

|Im ulrI < 3 X 10“® . (74) 


Eq. (73) gives at 90% C.L. 


This bound is at least two orders of magnitude below the bound obtained from /3 decay. 
Improving this bound in /3 decay requires a measurement of Dt < 10“®, which is an order 
of magnitude below the contribution of the ESI. 
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(a) Tree-level contribution to /3 decay 



(b) Loop contribution to the neutron EDM 


Figure 9 Example of scalar LQ exchange that contributes to (a) f5 decay at loop level, and to (b) 
the neutron EDM via an electroweak loop. The scalar LQ are denoted by R- and i?+, where ± 
refers to the weak isospin component (Ng and Tulin, 20121. 


The result above is obtained in a model-independent EFT approach, by introducing 
dimension-6 operators. The constraints apply to left-right symmetric models, exotic fermion 


models, and the R-parity violating MSSM (Ng and Tulin, 2012). Evasion of the bounds in 


Eq. (74) is only possible in either a strongly hne-tuned model or in a model in which the 


dimension-6 operators do not exist or do not contribute to either EDMs or /3 decay. An 
example of the latter are leptoquarks (LQs). LQs are particles with both baryon and lepton 
number, which can be either vector or scalar particles depending on their spin. These were 


previously considered “EDM-safe,” but in fact they are not (Ng and Tulin, 2012). LQs can 


contribute to jS decay at tree level, for example via the exchange of scalar LQs as depicted 


in Fig. 9a Leptoquarks also contribute to EDMs, but only the W exchange (Fig. 9b). Ng 


and Tulin (2012) show that these loop contributions are not suppressed by the light quark 


masses as was previously argued (Herczeg, 2001). Therefore, the constraints from 


EDMs in the LQ scenario are much more stringent than previously thought. 


Estimates of the limit on Dt in this scenario depend on the LQ mass and on whether 
light right-handed neutrinos exist. Assuming the existence of light right-handed neutrinos. 


Ng and Tulin (2012) found 


Im ulr = Dt/ao < 3 x 10 


-4 


'300 GeV 


(76) 
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while without them 


Im ttLR = Dt/ao < 7 X 10 


-5 


'300 GeV' 


(76) 


Ng and Tulin (2012) conservatively take miq = 300 GeV, which would give, assuming the 
existence of light right-handed neutrinos, Df < 3 x 10“^, a limit of the same order as the 
current /3-decay bounds. Nevertheless, improving the current (3 decay limit seems a difficult 
task, since there are many experiments ongoing or planned that aim to improve the bounds 


on the neutron EDM ( 

Altarev et al. 

2009 

2012; 

Baker et al. 

2011 

van der Grinten et al. 

2009 

Ito 

2007 

Serebrov et al. 

2009 

). In addition, strong bounds 

on the scalar LQ mass 


exist from the ATLAS and GMS experiments at the LHG. The bounds on their masses range 


from 607 GeV to 830 GeV, depending on the assumed LQ branching ratio (Olive et al ., 2014), 
which suggests much stronger bounds on Dt. 


2. Limits on R from EDM limits 


The R coefficient and the EDM are both P- and T-odd. EDM measurements in atoms 
and molecules limit both the electron EDM and BSM scalar and tensor electron-nucleon 


interactions. Khriplovich (1991) showed the relation between these electron-nucleon inter¬ 


actions and the electron-quark interaction of f3 decay. The scalar and tensor electron-nucleon 
interactions are defined by 


D = ^ ^ ICsNNei'^^e + CTNa^^Nei'^^a^^"'e 

N V2 


(77) 


where Cs {Ct) is the scalar (tensor) coupling and we have neglected pseudoscalar couplings. 


In Khriplovich (1991) and Khriplovich and Lamoreaux (1997) it was shown that the limits 


on Cs and Ct can be related to both Im and Im aiL) the couplings contributing to the 
R coefficient. 

The best current limit on nucleon scalar couplings is due to the EDM limit on molecular 


ThO, IG^I < 5.9 X 10 ® (90% G.L.) (Baron et al., 2014). The best bound on the nucleon 


tensor coupling, \Ct\ < 1.3 x 10 ® (90% G.L.), is derived from the EDM limit on atomic 


Hg (Ginges and Flambaum, 2004; Griffith et al, 2009). These couplings must be translated 
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Figure 10 Contribution of /3-decay coupling to the effective P- and T-odd electron-quark coupling 
through the exchange of the W boson. 


to quark couplings in order to compare them to the /3-decay couplings in Eq. 0- At the 
quark level, scalar and tensor couplings in the electron-quark (e-q) interaction are described 


by (Herczeg, 2003) 


Gp 

^ = ^ —Jp [ksq{ei'-f^eqq) -k kpqiei'y^aq)] , 

q=u,d V 2 


(78) 


where ksqikpq) is the scalar (tensor) coupling in the e-q interaction. The nucleon couplings 


can be translated into quark couplings by using the calculations in Herczeg (2003, 2005), 


which show that nucleon and quark couplings are of the same order of magnitude. Conser¬ 
vatively, we hnd that the ksq and kpq couplings are < 10“® (90% C.L.) 


Figure 10 shows that the electroweak corrections to the exotic /3-decay couplings con¬ 


tribute to the EDM e-u couplings, ksu and ksd- The effective P- and T-odd e-u interaction 


in Fig. 10 is estimated as (Khriplovich, 1991) 
—Gp cr 


a/ 2 47r 


In 


h 


Vud Im {2Al -k 24aL) ei'y^euu + u 


(79) 


where /i is the renormalization scale. Limits on the scalar electron-nucleon interaction Cg 
thus limit both Ag and ap. The effective e-d interaction only contains A^, and gives similar 
constraints. 


Comparing Eq. (78) and Eq. (79) we arrive at an expression for kgu and kpu- By us¬ 
ing kgu < 10“® (90% C.L.) and the conservative assumption that ln(/i^/m^) = 1 as in 


Khriplovich (1991), we estimate that at 90 % C.L. 


-5 


|Im %£,! < 10 
|Im < 10“® 


(80a) 

(80b) 
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Im aiR 

\m Al 

Im 

(3 decay 

4 X 10“^ 

6 X 10-2 

3 X 10-5 

EDM 

3 X 10“® 

10-5 

o 

1 


(3 X 10-^) 




Table VII Comparison between /S-decay limits on imaginary couplings and constraints from EDMs. 
The bound in parentheses is derived in a model with leptoquarks and right-handed neutrinos. For 
the /I-decay coefficients we use gs = 1.02(11) (Gonzalez-Alonso and Camalich 20141 and qt = 
1.047(61) (Bhattacharya et al. 2014) and the R-fit mentioned described in Sec. IV.B, Constraints 
are at 90% C.L. 


Both bounds are at least two orders of magnitude better than those obtained from the R 
coefficient in (3 decay. 


C. Conclusion 


Table VII summarizes the limits on imaginary couplings. Bounds obtained from EDMs 
are several orders of magnitude better than current bounds from T-violating /9-decay coef- 
hcients. The many ongoing efforts in the EDM held will strengthen the EDM bounds even 
further. 

The D coefficient should be measured with a precision of 10“® to improve the current 
EDM limits. Such a measurement is below the ESI interactions, and would require precise 
knowledge of the ESI for the used isotope. Measurements of the D coefficient are considered 


as part of a larger effort to measure 11 coefficients ( R) in neutron decay (Bodek et al. 


2011). Measurements of D are also considered in nuclear decays (Behr et al, 2014 Lienard 


2014). The E coefficient in Eq. (67) depends on the same BSM coupling as the D coefficient 


and is thus subject to the same EDM constraints. 

It might be possible that the connection between EDMs and (3 decay is diminished in a 
specihc new-physics model, when such a model is strongly hne-tuned. For the D coefficient 
examples are leptoquark models. Conservatively, this model relaxes the EDM constraint 


by maximally two orders of magnitude to \Dt\ < 3 x 10 (Ng and Tulin, 2012). This is 


of the same order as current [3 decay limits. Direct bounds on leptoquarks from the LHC 
experiments already suggest a stronger bound. Besides that, new bounds on the neutron 
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EDM are also expected before any new D measnrement could realistically be done. This 


would further improve the bounds in Table VII 


Improving the current bounds on Im ai requires a measurement of < 10 which is an 
improvement of the current result by more than 3 orders of magnitude. An R measurement 


in °Li is ongoing at the Mott polarimeter for T-vioIation (MTV) (Totsuka et al, 2014). 
Specihc models may again weaken the connection between jS decay and EDMs. Such models 


would have to be strongly hne-tuned. For example, Herczeg (2005) showed that in R-parity 


violating SUSY (Herczeg, 2005) such a cancellation would have to occur over 3 orders of 
magnitude. Such a severe cancellation is highly unnatural. Besides EDM limits there are 
also strong limits from the ratio = r(7r —)■ eu)/r{7i —)■ /ii/), which give Im < 4 x lO”"^ 


(Herczeg, 1995, 2005). 

Our EFT approach only applies when new physics can be parametrized by the heavy scale 
of new physics. If new particles are very light, the EFT approach does not apply anymore. 
However, the absence of experimental evidence for such light degrees of freedom supports 
the validity of the EFT approach. We therefore conclude that new measurements of the D 
and R coefficients should take the EDM bounds into account, and stress that the bounds 
can only be evaded in specihc and strongly hne-tuned models. 


VI. LORENTZ VIOLATION 

We will now review the new held of searches for the violation of Lorentz symmetry in 
the weak interaction. Recently, it was found that /? decay ohers unique possibilities to test 
Lorentz and/or CPT-invariance in the weak interaction, in both the gauge and the neutrino 
sector. We discuss these two sectors separately. 


A. Gauge sector 

In the gauge sector, Lorentz violation can be studied in a general theoretical framework. 


developed to study allowed and forbidden jS decay and orbital electron capture (Noordmans 


et al, 2013a||b Vos et al., 2015b). This framework considers a broad class of Lorentz- 
violating ehects on the W boson, by adding a general tensor to the Minkowski metric. 
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At low energies, this modifies the hh-boson propagator to 


^ , (81) 

where g'^'^ is the Minkowski metric and Mw is the hh-boson mass. Vertex corrections are 
described by 

-tT = -tg{g^’^ + Xn • (82) 


However, snch vertex modifications also requires the modification of the electron and neu¬ 


trino spinors (Noordmans et al ., 2013b). We restrict ourselves to propagator corrections, for 
which hermiticity of the Lagrangian implies that x]iu{p) — Xufi{—p)- In terms of the SME 


discussed in Sec. II.B one finds, at lowest order. 


A “ 2g^<pW ^ 


uQpQo 


'^w 1^2 


(83) 


where q is the momentum of the W boson and g is the SU{2) coupling constant. 


Bounds on y have been derived from allowed (Bodek et al. ^ 2014 Muller et al., 2013 


Wilschut et al ., 2013) and forbidden /? decay (Noordmans et al., 2013a), pion decay (Altschul 


2013 Noordmans and Vos, 2014), muon decay (Noordmans et al.. 2015), and nonleptonic 


kaon decay (Vos et al, 2014). Here we discuss allowed and forbidden jS decay. 


1. Allowed (3 decay 


For allowed {3 decay, Noordmans et al. (2013b) derived the Lorentz-violating differential 


decay rate using the modified W-boson propagator in Eq. (81). The complete expression 
is given in Eq. (IBII). Lorentz violation gives many additional correlations, since the ob¬ 


servables (momentum and spin) can now also couple to the tensor y. In /3 decay, a vari¬ 
ety of correlations can be used to access different (combinations of) y components. The 
necessary expressions can be derived by integrating over one or more kinematic variables. 
Momentum-dependent terms are always suppressed by some power of a heavy mass (M^ 
in the least-suppressed case), and can therefore be neglected given the current experimental 
precision. Neglecting momentum-dependent contributions to the propagator, the relation 
Xliuip) — Xvp.{~p) implies that y can only be real and symmetric or imaginary and anti- 
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symmetric, i.e. Xr^ = X^r i xf = Xi^ = 0) x’‘r = Xr^" = ~xf- The subscripts 

r and i denote the real and imaginary parts of y, respectively. This leaves 15 independent 
CPT-even components of that need to be measured. 

With this simplihcation and in the absence of tensor polarizations, the decay rate is 


Noordmans et al., 

2013b; 

Vos et al. 

2015a 


dW = 


F{±Z,E, 


— \Pe\Ee{Ee — Eq^ dEedfledflu^ |l + (2a — c)Xr^ + (—(2a — c )y(!^ + 2^y(^ 


E, 


pjpj 

EeE^ 


(a + c' + 2ay°°)5j7 - ^Xr] “ (2a - c')y' 
.1 


Os 


iPu X PeY 

EeEy 


^ i-2Lt + [(21 + Bxf)5u - 5y“] ) + 


JE, 


+ ^ ((—2a + c )y°-^ — 2gx^^ + {B + Axr^)5kj — Axr^ — Bx, 


{J^)pI r 




kj T^^.osiPux {j)y 


JE,, 


(84) 


where (J) is the expectation value of the spin of the parent nucleus, y^ = and 


Latin indices run over spatial directions. The last line of Eq. (84) contains only the neutrino 


momentum or the neutrino momentum and the nuclear polarization, and can therefore 
mostly be ignored. In fact, the neutrino correlations give access to a similar combination 
of y components as the electron correlations. The latter are considerably easier to obtain, 
and we will further only consider the electron correlation^ The coefficients a. A, and B 
are the standard /3-decay coefficients listed in Appendix the coefficient c' is a modified c 
coefficient. The coefficients with a breve (") multiply Lorentz-violating coefficients. They 
are given by 


c = 

9 
L 

d = 




1 + 


rA 


J'J 


+ ic 

Wot- , 


l+p2 ' 2 
l+p2 ’ 
1 + 


+ ■> 


(85) 


® In electron capture, the neutrino correlations play an important role (Vos et al. 2015bI. 
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where the upper (lower) sign refers to decay, and Xj'j and Aj/j = Aj/j 


are the standard /9-decay coefficients given in Eq. (A6) and Eq. (A7), respectively. The 


coefficient d vanishes for non-oriented nuclei and for nuclei with J' = J = ^. 

The effect of Lorentz violation in (3 decay can already be studied by measuring the 
dependence of the decay rate as a function of the direction of the emitted (3 particles. The 
modified Fermi decay rate integrated over neutrino energy and direction and summed over 
electron spin is 


dWp = dW^{l + 2x; - 2a' 




( 86 ) 


while for Gamow-Teller transitions of randomly-oriented nuclei 


dWoT = dtr" (i - + ?(x'“ + Xi)^ 


(87) 


where 


dW" = - EXF(±Z, E,:)dE, dni 


( 88 ) 


The component Xi can also be accessed by measuring the Gamow-Teller decays of polarized 
nuclei as a function of the spin direction, 


dWcT = dW^ 


{fy 


1 - 3 AT T Aj/jx' j 


(89) 


As an example of a mixed decay, one has for the neutron a = —0.11, A = —0.12,5 = 0.98, 
and g = L = = 0.27. Integrated over the neutrino directiorj® 


dW = dW^ll- 0.21x°° + (0.21x°' + 0.55x^)^ 


j 


-0.55x^ + (-0.12 + 0.98x“)|^ - 0.98xt"|^ - 0.12x»- ^ ^ ^ ^ (go) 


Ik Pe 


,0s (Pe ^ 


JE, 


Equation (84) depends on SM parameters, which are often not known better than at the 


1%-0.1% level. This dependence on SM coefficients can be avoided by measuring asymme¬ 
tries that do not depend on the accuracy of the SM coefficients. The Lorentz-violating part 


of Eq. (86) can, for example, be accessed by measuring the decay asymmetry of a Fermi 


This formula corrects Eq. (38) in 


Noordmans et al. 


(see also Appendix [Bj) . 
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transition with the /3 particles measured in opposite directions, 


W+-Wp 
W^ + Wp 


p: 


(91) 


where f3 = \pe\/Ee and Wp is the rate of (3 particles measured in the ± pe-direction. Simi¬ 
larly, the decay asymmetry in Gamow-Teller decays is 


Act — 


Wqt ^gt 

WSt + War 


i/?(x“ + x')pi 


(92) 


The coefficient x can also be obtained by measuring the spin asymmetry in a pure Gamow- 


Teller transition 


_ Wjr - W^r 
~ w^r + Wgt 




(93) 


't'f I 

where Wqp are the integrated decay rates independent of jS direction, but in the inverted 
polarization direction j, and P is the degree of nuclear polarization. A is the f3 asymmetry 
coefficient (for Gamow-Teller decays A = ^Xj'j). The remaining components of y require 
more complicated measurements that involve at least two observables. The decay asymmetry 
between the spin and the (3 particles can, for example, be measured from 


Ajjs — 


M-'lM-'r 


wlwi + Wi^wt 


'^'iWt ^ .2 PIHAx?< 


slk 


+ BXr)i‘p, 


(94) 


where Wl,r is obtained by measuring the (3 particles in the opposite pe direction, while the 
nuclei are polarized in the t (i) opposite j-direction. Similarly, can also be obtained 
by measuring the decay asymmetry between the neutrino and electron in perpendicular 
directions. 


The spatial directions of y are dehned in the laboratory frame and their absolute orienta¬ 
tion will depend on the orientation of Earth. It is therefore necessary to choose a standard 
absolute reference frame, for which the Sun-centered inertial reference frame is commonly 


chosen (Kostelecky and Russell, 2011). The movement of this reference frame can safely be 
ignored. The transformation of in the laboratory frame to the Sun-centered frame, in 


which we denote y^*^ by is (Noordmans et al, 2013b) 


X 




(95) 
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Sidereal days 


Figure 11 (Color online) Illustration of the oscillation of the asymmetry Ap in Eq. (91), for 
= 0.1 and C = 45°. Three different detections directions of the /3 particles are depicted. When 
13 particles are detected parallel (||) to Earth’s rotation axis there is no sidereal variation (red line). 
The blue line shows the asymmetry when the jS particles are detected in the east-west direction (T) 
direction and black line when they are detected perpendicular to the Earth’s surface (ti)- Both 
show a sidereal variation, the latter with a constant offset. 


The transformation matrix is 




(96) 


0 cos C cos fit cos C sin fit — sin ( 


0 — sin fit cos fit 


^ 0 sin C cos fit sin ( sin fit cos ( j 


where C is the colatitude of the experiment and hi is Earth’s sidereal rotation frequency. In 
the laboratory frame, x points in the north to south direction, y points west to east, and z 
is perpendicular to Earth’s surface. The coefficients x\ can be transformed to 

and Xl, respectively. This transformation shows that the asymmetries Ap, Act, and Aj 
can oscillate with the rotational frequency of Earth. These sidereal variations of the signal 
are a unique signature of Lorentz violation, and can therefore be separated from other limits 
on BSM physics. A generic example of how sidereal oscillations can be observed is shown in 
for = 0.1. This example also shows that if the (3 particles are detected parallel 


Fig. 


11 


(II) to Earth’s rotation axis, the asymmetry will have no sidereal dependence (red line). The 
blue line shows the case where the (3 particles are detected in the east-west (T)direction. It 
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- 0.1 
- 0.2 

Figure 12 (Color online) Illustration of the possible sidereal variations of tensor Lorentz violation 
parametrized as x^rJ^Pei with = 0.1. The red line (line 1) shows the modulations when j is in 
the z (up-down) direction and Pe in the y (east-west) direction. The blue line (2) is for j in the z 
direction and the (3 particles detected parallel to Earth’s rotation axis. The black line (3) shows 
the modulations when both j and pe are in the east-west direction. 

has no offset because it is measured perpendicular to Earth’s rotation axis. The black line 
gives the asymmetry for (3 particles detected in the up-down (t|) direction perpendicular 
to the Earth’s surface {z direction in the labframe). It shows a sidereal oscillation on a 
constant offset. Detection of the (3 particles perpendicular to the rotation axis is preferred, 
since an offset could be the result of systematic errors in the measurement. 

Tensor contributions involving lead to terms that may oscillate with twice Earth’s 
rotational frequency. Figure illustrates three possible scenarios for an asymmetry that 
depends on The red line shows the modulations when the polarization is in the 

up-down direction, while the (3 particles are detected in the east-west direction. The blue 
line shows the modulations in the same polarization direction, but when the (3 particles are 
detected parallel to Earth’s rotation axis. It shows an oscillation with the period of the 
sidereal rotational frequency on top of a constant offset. The black line shows an oscillation 
with twice the period of the sidereal frequency. It arises when both the polarization and the 
(3 particles are detected in the east-west direction. 

In allowed f3 decay, Lorentz violation was for the hrst time tested in polarized ^‘^Na 
, by measuring the spin asymmetry Aj (Eq. ([9^). ^°Na hrst decays 
with a /d'*' —J- 2^ Gamow-Teller transition, followed by a 7 decay of the daughter nucleus. 

The parity-odd (3 decay was used to determine the polarization P by measuring the (3 


(Muller et al. ^ 2013 



0 1 2 


Sidereal days 
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Figure 13 Schematic setup of the rotating experiment 


iNewiiian anu vviesnei 




asymmetry (Muller et al, 2013). The parity-even 7 decay was used to measure the lifetime 
7 -t(l) determine the 7 asymmetry 


rt 


yl.y — 


7-t -|- rl 


= P^Xi ■ j , 


(97) 


where the polarization direction is in the j direction. To reduce systematic errors, the 
polarization direction is preferably in the y (east-west) direction. The analysis of the setup 
in this direction places bounds of order 0(10“^) 


(Sytema et al ., 2015 


Lorentz violation has also been searched for in polarized neutron decay (Bodek et al. 


2014). Two different asymmetries, that depend on the nuclear polarization and the {3 direc¬ 


tion, were measured and are currently being analyzed. The asymmetries depend on combi¬ 
nations of Xi and Xr and preliminary bounds are 0(10“^) (Bodek et al.. 2014). This setup 


probably also allows for a measurement of Ajy dehned in Eq. (94). Such a measurement 
would measure the so-far unconstrained coefficients x?- 
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2. Forbidden jS decay 


“Forbidden” (slow) transitions are suppressed with respect to allowed transitions, be¬ 
cause the lepton pair carries away angular momentum. Theoretically, the simplest of these 
transitions are the unique hrst-forbidden transitions (AJ = 2), since they depend on only 
one nuclear matrix element. Because Lorentz violation includes rotational violation, it also 
implies the violation of angular-momentum conservation. Forbidden jS decays are then more 
sensitive to rotational invariance violation in the weak interaction. In the 1970s, two ex¬ 


periments were performed with this motivation. Newman and Wiesner (1976) searched for 


anisotropies in the angular distribution of /3 particles in hrst-forbidden decay. Ullman 
(1978) searched for sidereal modulations of the count rates for hrst-forbidden ^^^Cs (3 decay 


and second-forbidden ®®Tc (3 decay. The strongest bounds were found in the experiment 


by Newman and Wiesner (1976). In this experiment the /3-decay distribution of from 


a high-intensity source was measured in a rotating setup. Schematically, the setup is de¬ 


picted in Fig. 13 The rotation of the setup allowed for the determination of three decay 
asymmetries 

, ^Wr,-Ws , ^IY^-IYh. 


and 


Sou — 2 


Wn + Ws-We- IVw 


(99) 


Wn + Ws + We + Ww ’ 
where A, S', E, W are north, south, east, and west and W is the decay rate measured by the 
f3 particles in that direction. These asymmetries were htted with 


6 = Qq + ai sin(f 2 t -|- ^i) -7 02 sin( 2 f 2 t -|- 02 ) , 


( 100 ) 


to search for a sidereal time dependence and to reduce systematic errors. The extracted 


bounds on ao,ai, and 02 are (9(10 ®) (Newman and Wiesner, 1976). Noordmans et al 


(2013a) reinterpreted the data from Newman and Wiesner (1976) and Ullman ( jl978 ), after 
extending the allowed /3-decay framework to include higher-order terms in the multipole 
expansion, i.e. all possible forbidden decays. The modihed IF-boson propagator gives an 
unconventional contraction of the nucleon and lepton currents, such that angular momentum 
is no longer conserved. In the Lorentz-symmetric case, rotational invariance implies that 


63 

























AJ AvliGiTG /\ . 1 IS IjItg spin ch^jn^G of tho nuclons Bind t^iep is tlio totnl ctn^ulcir nionion~ 

turn of the leptons. In contrast, when contracting with transitions with AJ = Jiep + 1 
are possible, and when contracting with cilso A J = Jiep + 2 transitions are allowed. It 
is thus possible to have transitions in which the leptons carry away less angular momen¬ 
tum than in the Lorentz-symmetric case. Because the suppression of the forbidden decays 
is proportional to the angular momentum of the leptons, the Lorentz-violating terms are 
enhanced compared to the Lorentz-symmetric case. 


For unique hrst-forbidden transitions (Noordmans et ai, 2013a) 
dW 


dVtfdEf 


2 I 2 I 2 


aZ - 

PeR L 


ila (xrW - - \x\p' + x'rV 


( 101 ) 


where aZ/peR — Equation (101) shows that the Lorentz-violating contributions 


are enhanced. Higher-order forbidden decays do not have additional enhancement compared 


to the simpler first-forbidden transitions. Noordmans et al. (2013a) translated the bounds 


from Newman and Wiesner ( 1976| using Eq. (101). This led to strong limits on several 
combinations of Assuming no cancellations between coefficients, this results in the 


limits (Noordmans et al. ^ 2013a) 
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( 102 ) 


These are the strongest constraints on The only coefficients not constrained by for¬ 
bidden decays are yf. These coefficients can be studied in allowed fd decay by considering 


Eq. (94) or equivalent correlations. The bounds on y were also translated into bounds on 


the SME parameters (Noordmans et al., 2013a), providing strong direct bounds on the SME 


parameters and /c^vv dehned in Eq. (16) 
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3. Conclusion and outlook 


We have discussed the efforts to search for Lorentz violation in the weak interaction in 
forbidden and allowed (3 decay. The bounds from forbidden (3 decay are several orders of 
magnitude stronger than the current bounds in allowed (3 decay, due to the intense sources 


that were used (Newman and Wiesner 


Ullman, 1978). In allowed {3 decay, Lorentz- 


violating effects are not enhanced and matching the statistical precision of forbidden (3- 
decay experiments would require long-running experiments with high-intensity sources. An 
interesting alternative lies in orbital electron capture, where it is possible to use such high- 


intensity sources (Vos et al. ^ 2015b). 


Allowed (3 decay offers various correlations in which Lorentz violation could be probed. 
Observables can be chosen such that they give direct constraints on y compared to the com¬ 
bination of coefficients constrained by forbidden decays. Two relatively simple experiments 
that probe the (3 asymmetry in Fermi and Gamow-Teller decays, Ap and Aqt respectively, 
would give direct bounds on and Xi- These asymmetries could be studied parallel to the 


efforts to measure the /S-spectrum shape discussed in Sec. IV.D (Vos et al ., 2015a). Another 
interesting possibility is to exploit the 7 ^ enhancement of decay asymmetries by considering 


fast-moving nuclei (Altschuh 2013| Vos et al ., [2014, 2015a). The total decay rate in the rest 


frame of the nucleus is proportional to Xr Eq. ( 86 ) and Eq. (87)). For a fast-moving 


nucleus, the expression can be related to the Sun-centered frame with a boost. If the nucleus 
is moving ultra-relativistically in the v direction. 


= Jr ( a '7 + , 


(103) 


where 7 ^ is the Lorentz-boost factor and T, L, and K are coordinates in the Sun-centered 
reference frame. This relation was, for example, used to extract bounds of 0(10“"^) from 


pion decay (Altschuh 2013). For allowed (3 decay, /3-beam facilities, currently considered for 


producing neutrino beams ( ]Lindroos and Mezzetto , 2010), could be exploited. 

So far the coefficients xf remain unconstrained. In Fermi decays, this coefficient could 
be measured by considering the correlation XiKPe ^ PuY- The coefficients can also be con¬ 


strained by measuring the polarized /3 asymmetry Ajjs in Eq. (94). Such an asymmetry could 


probably be explored in the neutron-decay measurement pursued by Bodek et al. (2014) 
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B. Neutrino sector 


A different possibility to study Lorentz violation in (3 decay lies in the neutrino sector of 


the SME (Kostelecky and Mewes, 2004, 2012). Most interesting for (3 decay are the modihed 


versions of and dehned in Eq. (15). 

Unlike the gauge sector, the neutrino sector has been studied extensively in several ex¬ 
periments. Strong bounds exist from neutrino oscillations and time-of-flight measurements 


(Kostelecky and Russell, 2011). However, there are four operators that do not show up in 
oscillations and have no effect on the neutrino group velocity. These operators are called 


“countershaded” (Kostelecky and Tasson, 2009). Recently, Diaz et al. (2013) showed that (3 
decay has a unique sensitivity to these operators. The four countershaded coefficients are 
denoted by The operators are dimension 3 and CPT-odd. These coefficients modify 
the neutrino dispersion relation and the available phase space of the neutrino, which affects 
f3 decay in two ways, in the (3 end-point and in the correlation coefficients. 


1. End-point in f3 decay 

The /3-spectrum end-point is very sensitive to the neutrino phase space and to the neutrino 


mass (see also Sec. IV.C). Independent of the neutrino mass, the countershaded neutrino 


coefficients also shift the end-point, as can be seen from the modihed decay rate (Diaz, 2014 

dW 


Diaz et al ., 2013) 


dT 


[AT + iTw)" - kml , 


(104) 


where AT = Tq — T^, = Ee — me is the electron kinetic energy, and Tq is the end-point 

energy for = 0. 5Tlv is the Lorentz-violating modihcation, which depends on sidereal 
time. Independent of the neutrino mass, a bound on the countershaded coefficients can be 


set by using the available data of the Troitsk (Kraus et al.. 2005) and Mainz (Aseev et al. 


2011) experiments, see Diaz et al. (2013). Since these experiments collected data over a long 


period of time, all the oscillations average out and only the time-averaged Lorentz-violating 
coefficients can be constrained. Therefore, only two of the four countershaded coefficients 


could be bounded. Conservatively, the analysis of Diaz et al. (2013) gives bounds of order 
0(10“®) GeV. These limits improve and complement previous limits. A dedicated analysis 


of the data of the Troitsk, Mainz, or the expected KATRIN (Otten and Weinheimer, 2008) 
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experiments could improve these results. If the data analysis also takes into account the 
sidereal time, bounds on all the countershaded coefficients could be set. 


2. Correlation coefficients 


The Lorentz-violating neutrino coefficients of Eq. (15) also modify the neutrino spinor 
solutions. Near the end point, this modihcation can be neglected because the phase space 
dominates. However, the derivation of the complete modihed decay rate requires both the 
modihed spinors and the phase-space modihcation. The modihed neutrino phase space is 
~ (E‘^ — 2Euaf^)dEydrtv. The modihcation of the spinors requires the replacement of 
Piy by py = {pu -|- — djjfpy), where is the isotropic component. The modihed neutron 

decay rate is 


dW 
dVtp dVt,, dT 


~ E{Z, E)\ffi\Ee{El + 2EffiTuj) 


(105) 


The neutrino coefficients modify the decay rate in a similar way as y does, since there are 
now additional correlations between J and and 

The countershaded coefficients could, for example, ahect the (d-v correlation. The jd-v 
correlation can be measured as an asymmetry, dehned by 


_ N+-N_ 
~ N+ + nS 


(106) 


where A^+(iV_) is the number of decays in which the neutrino and electron are emitted 
(anti)parallel. The Lorentz-violating neutrino coefficients modify this correlation coefficient 


to ( [Diaz 2014) 


a = a\m + ./I + 

V tt E,. 


^of no 1 


(107) 


where the coefficients should be transformed to the Sun-centered frame and would depend 
on the sidereal frequency of Earth. 


No experiment has searched for these variations, but Diaz et al. (2013) estimate that a 
0.1% measurement of a would limit the countershaded coefficients at the level of 10“® GeV. 
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Similar, for a 0.1% measurement of the correlation coefficient B, the limits are estimated at 
0(10“®) GeV. A dedicated experiment measuring either a or B would thus provide inter¬ 
esting new bounds on Lorentz-violating parameters in the neutrino sector. Note that x and 
have a similar influence on the decay rate. In a dedicated experiment both coefficients 
might influence the asymmetry. A measurement of Eq. (106) might also be sensitive to xf 
and Xiy depending on the experimental setup. 


C. Conclusion 


To summarize, /3 decay offers a unique way to study Lorentz violation in both the gauge 
and neutrino sector. The large variety of correlations allows for direct measurements of 
different components of y, while in the neutrino sector (3 decay allows for the study of 
countershaded coefficients. 


In the gauge sector, strong bounds on the order of 10“® - 10“^ exist from forbidden (3- 
decay experiments. Unconstrained are the coefficients x?, which can be accessed in f3 decay 
by considering the interaction of x with two observables (Eq. ( [9^ . Improving the existing 
bounds requires high statistics and precise knowledge of the systematic uncertainties. Ben¬ 
eficial for this would be to exploit the Xr enhancement of boosted (3 decay or to consider 
electron capture. The real and imaginary part of x can be constrained by measuring the 


asymmetries in Eq. (91) and Eq. (92), respectively. Such an effort could be combined with 


measurements of the Fierz-interference term. 


Further, we have discussed the possibilities to improve constraints on Lorentz violation 
in the countershaded neutrino sector. In that sector no dedicated experiment has been 
preformed so far, but using available data from tritium gives bounds of order 10“® GeV. 
The parameters not constrained so far could be bound in /3-decay correlation experiments. 
Lorentz violation gives a unique signal compared to other BSM physics when searched for in 
a dedicated experiment. Estimates for 0.1% measurements of the coefficients a and B gives 
a constraint on Lorentz violation of 10“^ GeV, which shows the potential for these future 
experiments. 





VII. SUMMARY AND DISCUSSION 


In this review we addressed the current status and role of nuclear and neutron f3 decay in 
the search for physics beyond the SM. In these searches, the statistical precision is becoming 
increasingly important. However, systematic errors, despite improved detection methods, 
and higher-order corrections such as FSI, still appear to be the main limits. In the meantime, 
thanks to the evolution of EFT methods, constraints obtained in other helds weigh in, 
establishing bounds on the scalar and tensor contributions. This is illustrated in Fig. 
and Fig. where measurements at LHC (St 


5ec. 


IV.B) and limits from the neutrino mass 


(Sec. IV. CD give constraints that outperformed the /3-correlation measurements in the right¬ 


handed sector. This is quantihed in Table IV 


The study of fundamental aspects of /3 decay will be most fruitful in the study of left¬ 


handed scalar (Sec. IV.A.1) and tensor currents (Sec. IV.A.2), as these appear linearly 
in most observables via the Fierz-interference term. Fortunately, these interactions can 


be studied in parallel to precision studies of SM parameters (Sec. HI.B). For example, 
extracting the CKM matrix element Vud from superallowed Fermi transitions has, as a by¬ 
product, the most strict limit on left-handed scalar interactions. Lacking still is a similar 
limit on tensor contributions. An interesting option to obtain such a bound could come 
from measuring the detailed shape of the (3 spectrum in Gamow-Teller transitions. Also the 
potential of mirror transitions, both for obtaining tensor limits and for obtaining a value 
for Vud independent of the superallowed Fermi transitions, has been recognized. In Table |V] 
we indicate the precision required to impose new bounds on left- and right-handed scalar 
and tensor currents. Measuring the Fierz interference term in /3 decay remains competitive 
in determining bounds on left-handed coupling constants. In contrast. Table V shows that, 
for right-handed couplings, the limits from LHC and the limits derived from the neutrino 
mass are by far superior to the best bounds derived from the /3i/-correlation a, and future 
experiments in f3 decay are unlikely to reach this precision. 

Concerning the most fundamental measurement of T-violation, we discussed in Sec. |V] 
the strong bounds on the triple-correlation coefficients D and R derived from the limits on 
permanent EDMs. These bounds are summarized in Table VH[ Not only are the bounds 
from EDMs several orders of magnitude stronger than those of f3 decay, but the EDM limits 
also have a large potential to improve faster than those from /3 decay. One reason is that 
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EDMs can be measured in stable or long-lived particles, but also because of the widely 
perceived urgency for improved limits in this sector. 


A new twist to the discussion of symmetry violations in jS decay has been added, since 
(3 decay also offers an interesting sensitivity to Lorentz violation in the weak interaction. 
In Sec. |VI[ we reviewed these limits for the hrst time. Because the discrete symmetries C, 
P, and T are each violated in the weak interaction, this interaction is a promising portal 
to search for new physics when considering CPT violation and thus Lorentz violation. The 
familiar /3-decay correlations are now extended to include correlations between spin and 
momentum and a Lorentz-violating background tensor. Consequently, spin and momentum 
will appear to have preferred directions in absolute space, resulting in unique signals that 
can be distinguished from other BSM searches. 

In weak decays the LIV has been parametrized with the complex tensor y. The bound on 


most components of this tensor are of the order of 10“° to 10“® (Sec. VI.A.2). Fine tuning 
between the tensor components allow to weaken these bounds. Relatively simple new exper¬ 
iments can improve these bounds using very strong sources, also removing the possibility of 
hne tuning. Obtaining sufficient high counting statistics is the main challenge. The searches 
for Lorentz violation can be expanded in a parallel effort with the more traditional searches. 
Alternatively, one can study (3 decay in flight, exploiting the 7 ^ enhancement. In this respect 
there may be as yet unexplored possibilities related to semileptonic decays in high-energy 
physics. Because this held of research is relatively unexplored, both experimentally and 
theoretically, the best approach may still emerge. 

Improvements in theory and experimental techniques, as well as new radioactive-beam 
facilities, provide new possibilities to study fundamental aspects of (3 decay, both in the 
search for exotic interactions and in the search for Lorentz violation. These studies should 
be done by considering also the other searches in high-energy physics and precision physics 
at low energies. Nuclear and neutron f3 decay will remain an important topic on the research 
agenda. 
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Appendix A: Decay coefficients 


Yang 


Our formalism can be linked to the original work of Jackson et al. (1957a) and Lee and 
, wher^ 


fJeff = pn e{Cs - 0^75)1/6 

67^(Oy - Cy'y^)Ue 
+ Pl^l^n 67;,( 0(1 - CAlb)Ve 

+ ]^pa^''n ea^uiCr - 0^75)1/6 . (Al) 

This notation can be related to our couplings in Eq. by using the normalized couplings 

Q = A , (A2) 

and 

Cv = gv{(^L + cir) , 

C'v = gv{o.L — CIr) , 

Ca = ~\gA\{.0,L + 5 

C'a = ~\gA\{,af — a^) , 

Cs = gsi^L + Ar) , 

C's = gsi^L - Ah) , 

Ct = ‘^grio'L + Oh) , 

Cf = 2gT{aL — Q!h) • (A3) 

^ Using our definition of 75 and neglecting pseudoscalar couplings. 
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For simplicity we have defined 


OL = O'LL + O'LR , 

= O'RR + <lRL , 

<^'l = ~ ®- L-R 5 

/ _ 

O-R = CtRR — RrL , 
= ^LH + ^LL , 
+ ^RR ■ 


(A4) 


We write ai and or as in Eq. ( 0 . because 75 = . The coefficients ae 5 ,A^s, 


and tte are related to the e coefficients in Cirigliano et al. (2013b) and Naviliat-Cuncic and 


Gonzalez-Alonso (2013) by using 


{O'LLi O'ER, GIRLi O^RR, All + Alr, Arr + Arl, aL, Or} — {1 + Gl, Gr, Gl, Gr, Gs, Gs, 2gt, 2gt} ■ 

(AS) 

A full list of correlation coefficients in allowed jS decay including Coulomb corrections is given 
in Jackson et al. ( 1957a[b ). Here we will give the most important decay coefficients in terms 
of couplings defined in Eq. 0 . We emphasize that only b, B, and N depend linearly on 
scalar and tensor couplings. We define X = \gA\/gv > ^ and neglect Coulomb interactions. 
The spin factors are 

/ 

J ^ J' = J - 1 

J ^ J' = J , (A6) 

J ^ J' = J + 1 


Xjij — 


1 , 


1 

j+i’ 


-j 
j+i’ 


and 


Arj = -(2J-1) 


J ^ J' = J - 1 
J ^ J' = J 
J ^ J' = J + 1 


(A7) 


(J+l)(2J+3)’ 

where J and J' are the spin of the initial and final nucleus, respectively. In the following 
equations, Mp and Mgr are the Fermi and Camow-Teller matrix elements, the upper (lower) 
sign refers to (3~ decay, and 7 = \/l — a‘^Z‘^, with Z the atomic number of the daughter 
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nucleus and a the fine-structure constant. 


Because of our normalization of the couplings in Eq. 0 we define 




2 


with 


^ ||aL|^ + IctR^ + p- (I^lT + I 

+ 2(y'yA^|MGr|^ |l®Lr + + 4p (ictiT + | 


(A8) 


(A9) 


Neglecting Coulomb interactions, the decay coefficients are (Jackson et al. 1957a[b) 


ae = 2gl\MF\^ \aL? + ^ + |kl«| 

I 9v 

+ + 4?t ^ |„^|2 

-J I 9 A 


(AlO) 


6? = ± 2g^7 2|Mp| 


|2^ 

9v 


Re{ALal) + Re(Aija^) 


\9a\ 


Re(aLa^) -|- Re(aRa(^) 


(All) 


ce = 2glX^Aj,j\MGT\^ -K\^ - la'^r + 4^ + l«i?P 

I 9a 


(A12) 


Ae = 


± 25 (yA"‘|MGT|"‘Aj'j 4^ \aL\^ — IttR^ 

9a 


+ ‘2gy\6jij\MF\\MGT 
+ 2 


J {^ 9t9s 


<4 + 11 \9a\9v 

I^LlI^ — loLi?!^ — loHiJp + \0‘Rl\‘^ 


I M 2 I / |2 

PlI +ijl 


Re(Aia2) — Re(Aj:ja 


R) 


(A13) 
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— 2(7^/A \Mgt\ Aj/j 


-Aqt rUeJ 
\9a\ Ee 


Re{aLal) — Re(aija'^) 


± 


9a 


\(Al\ 


l«ijl 


± 


J 12 


4P 


25'y A 5 j' j I Mp 11 Mgt 


All 


\J ^ 

^9t9s 

1 J + l\ 

^9v\9a\ 

1 |2 1 
l^iRRl + 

1 |2 
\CiRL\ 


Re(R2,Q;2) — Re(Rija 


R) 


± 


rrie'y 


‘^9s 

9v 


o>lr\ 

Re(ALa'^) - Re{ARap) 


+ 


‘^9t 
\9a\ ■ 


Re^aial) - Re(aija;^) 


, (A14) 


— 2gYX6j'j\Mp\\MGT\ 


J I ‘^9t9s 
J + 1 Ifi'vIfi'Al - 


\m{ALa*p) + Im(ARa;^) 


+ 2 


Im(aLa'£) + Im(aija'^) 


(A15) 


Ri — E 2gyX \Mgt\ Aj'j- 


~4^t 

\9a\ 


‘^9v^^J'j\Mf\\Mgt 
A gr 


\9a\ ■ 

and 

= 2fi'yA^|MGT|^Aj/j 


J + 1 

Im^apal) - Im^apap) 


Im^apap) - Im(a;^a'^) 

Im^Apal) - Im(Aija'^) 


J 


2 X / mel 


1 E, 


Wl\‘^ + Wr\‘^ + 4^ 

9a 


I |2 I I |2 

k^L + k^_R 


' ^2 


Re(aLa2) + Re(aj?a^) 


t/A L 

‘^9v^^J'j\^f\\Mgt\ 

Agp 


J 


J + 1 


\9a\ - 


Re{aLCi*i) + Re(aija^) 


± 


I 12 I 12 I 12 I 

IOLlI — \0‘Lr\ + IOHrI — I^RL 


‘^iReiApa'D + ReiARa%) 
9v I 

'yrue / AgT9s 
Ee \9v\9a\ 

2 


Re{Apa*p) + Re(Aj:ja 


R) 


(A16) 


(A17) 


The longitudinal electron polarization is (Jackson et al, 1957a|[b) 


P = 




l + 6(g) • 


(A18) 
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with 


Ge= ±2\MF\^gl[^\\AL\^ 
\.9v 


-|ai|2 + |a^p 


± 2\Mgt\^ 


f ^9 t 
I 9a 


2 r 


ttL 


|«k| 


12 


7 ' |2 


(A19) 


The neutron lifetime (Eq. (40)) depends on Vud, which is extracted from the O’*" —)■ O’*" 
superallowed Fermi decays. However, the extracted value of Vud might also depend on new 
physics. Taking into account this possibility, 


t„. = K 




9v^ ^ ^ ^ 


1 + + 1^1 + 3V(1 + 4|ai + 4|a|) + (2^A, - ' 

(A20) 


where {nie/Ee)^ is the inverse average energy of the superallowed decays, calculated by 


Pattie et al. (2013). The constant K is (Pattie et ai, 2013) 


K = 


271^ 


miUl + AHc)G|Ci 


= (1.9342 ±0.002) ■ 10 


-4 


(A21) 


where /„ = 1.6887(2) is the statistical rate function (Towner and Hardy, 2010) and Arc are 


the SM electroweak corrections (Czarnecki et al., 2004). 


The SM expressions can be obtained by setting qll = ^ and neglecting all other couplings. 
Dehning p = \gA\MGT/9vMF, the remaining SM expressions are 


1 - pV3 

i±p^ ’ 

(A22a) 

TAj'jp^ + 2(5j/jy^J/(J + l)p 

Asm - 1,2 ’ 

1 + 

(A22b) 

EXjijP^ + 25 jij^ J/{J+ l)p 

EsM — 1,2 ’ 

1 + 

(A22c) 

Gsm = t1 , 

(A22d) 


while all other coefficients vanish. For neutron decay, p = \/^\gA\ and J = J' = \. 
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1. Linear terms in B 


The B coefficients contains terms linear in exotic couplings. Neglecting quadratic cou¬ 
plings, we can write 


Bi = ±Xj,jp^ + 2p5j,j^J/{J + l) + hsi , (A23) 

where 


bsi = 


2 X ^9T 

— P Aj'j-j—j-Ke ul 


\gA\ 


T Sj'jP\ 


J 


J + l 


2 fi'ST 5 /I I 

-Ke Al + 1 —rK.e ai 


9v 


\9a\' 


(A24) 


Most B measurements measure 


B = 


Bsm + 


1 - 1-6 


(A25) 


For pure Gamow-Teller transitions, with p —)■ cx), Bqt = ^^j'j and the linear dependence 
cancels. For neutron decay and assuming real couplings, Bsm — 1 and 


~ / 7Tl 

B — Bsm + — ^Bsm) 


2(A +A2 

_ 

1 + 3A2 

~ 1 -|- 


me7\ —A — 2A^ -I- 3A^ 
~^/ (1 + 3A2)2 

rrie'y 


r,9s . , A 9t 

2 — Al + 4:-—-as 

. 9v \9 a\ 


n A 9s A , n n 9t 

0.1— Al -h 0.2^ —-as 


9v 


\9a\ 




(A26) 


For comparison, for neutron decay, the Fierz interference term is 


b 


neutron — 


2^Al - 

9v ^ _ \9a\ ^ 

1 + 3A2 


~ 0.35—Ai - 3 . 3 /^aL . 
9v \9a\ 


(A27) 
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For the measured A coefficient in neutron decay, with AgM ~ —0.11, 

^neutron = ^SM =F ^SM(0.35 — . (A28) 

gv \gA\ 

So for neutron decay, B actually has a reduced sensitivity to scalar and in particular tensor 
terms compared to for example A. 


Appendix B: Lorentz violation 


The Lorentz-violating f3 decay rate including Coulomb corrections and electron spin, to 


hrst order in is (Noordmans et al. ^ 2013b) 


dW = 


{271 


-EePe{Eo - EefE{±Z, EA^dE^dAedA^ 


X j I It 
+T 2 


JE, 


Pe • 5 e 
Ee 

+ 


+ ^ jE^ 


E,, J 




E„ 


{Ee - ^nie) {Pe ' <?e) \ pi 7"ie 


E^ - ml 


I / / -1 V— ^ / \£' ^ c- / 1 /-'e / ' ' ^6 /V / /-I 9/^ '^ \ I 

+ I I 1 T ^ ^ 1 -4- T —< T ^Vl - lAPe X aeY 


X 


_ ZjP’P-jp , , J^PI , 




2^^ j 2- p ^ ■ 2 V- ' ->E' ' E, ' J 






JE,, 


+ 


E,. 


(Bl) 
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where 7 = Vl — The Lorentz-violating constants ar# 


? 

Wi = 


t = {a- ic')xr , 

-xx°r + Hxf -Xi) , W2= k{xf -xT)- LXi , 


= — 


'^Xr^ + + Xi) ) 


rjikm 3 J km 

-i 1 2 ^ X.r 


_ 

-‘-3 — 

ilk 


^ v/v.r /V/- / /v-i / o 

= \Ax^^5^^ + T(xf + xf + Xfe^^^) , 

^ 2 : + - {xxT + 9xfy^^ - 9{Xr + Xr ) , 

Tf = - L{kr - xfe^^') - k{Xr - xfy^) , 

kr = L (xf 5^™ - xfe^^^ - xf + xT^ - xfV”^^) 

+;x - xf5^^ - x^/y + (xr + xr 

/ /y^ Jyy /I O / / (V^ O J J'y rt /W-I 7^ t^J ^ /w^ / t^l/y /^ (V^ /I t^Ly] 


ms ^sjl 


XT 


_ 1^/ ^ml^kj ^ 


(B2) 


where r and i denote the real and imaginary part of x^^, respectively, x’' = ^i^k^mk^ ^z 
denotes the electron momentnm in the /-direction. a,A,B, and ^ are the standard /9-decay 


coefficients, given in Eq. (A22), the other coefficients are 


1 ~P 

^ 1 I 2 ’ P 1 I 2 ’ 

1 + p"' 1 + p"' 


9 = Kl -a:)(l + |AjjO , 


c' = (1 - x)Ajj> , 

k = - yJy ^6. 


J + 1 


jj' 


L = ±1- 




1 + p^ 


rA 


JJ' 


(B3) 


where npper(lower) signs refer to {3 (/d"*") decay. The coefficient Aj/j is given in Eq. (A6) 

A _ A {{J-3?)-\j{J+l) 

Aj/j = Aj/j^- TT7r-AT\ - 


and Aj/j = Aj/j^- j( 2 jk) -’ d^ - 
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